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AERL 


ABSTRACT 


The  possibility  of  film-cooling  the  B-wall  (insulator  wall) 
of  an  MHD  generator  duct  has  been  investigated  theoretical ly .  A 
potassium  seeded  nitrogen  plasma  as  the  primary  and  secondary  fluids 
is  considered  and  the  compressible,  turbulent  boundary-layer  equations 
are  solved  numerically  using  the  implicit  finite-difference  scheme  of 
Patankar  and  Spalding.  This  numerical  method  has  been  modified  in  a 
general  way  to  include  the  effects  of  electromagnetic  body  forces  such 
that  it  can  be  applied  to  other  than  film-cooling  flow  problems  which 
involve  the  action  of  such  body  forces.  The  equations  solved  are 
overall  continuity,  momentum  and  energy  as  well  as  electron  conserva¬ 
tion  and  energy  to  account  for  electron  thermal  non-equilibrium  and 
finite-rate  ionization  in  non-equilibrium  generators.  A  mixing-length 
scheme  to  model  MHD  damping  of  turbulent  fluctuations  is  proposed  by 
exploiting  the  analogy  which  can  be  drawn  between  MHD  and  wall  suction 
boundary  layer  flows. 

Results  show  that  the  adiabatic  film-cooling  effectiveness 
is  significantly  increased  and  that  MHD  turbulent  damping  noticeably 
increases  the  effectiveness  even  further  but  to  a  lesser  extent. 

An  optimum  magnetic  field  intensity  exists  for  which  maximum 
film-cooling  effectiveness  is  obtained  for  a  particular  flow 
condition.  The  MHD  film-cooling  process  is  found  to  be  as  sensitive 
or  less  to  initial  flow  parameters  and  geometry  as  has  been  observed  in 
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ordinary  hydrodynamic  flows,  but  initial  magnetic  field  perscrip- 
tions  must  be  as  true  to  physical  reality  as  possible.  Non-equilibrium 
effects  are  found  to  produce  near-wall  current  density  increases  which 
influence  the  film-cooling  effectiveness  in  a  complicated  way,  the 
net  outcome  of  which  results  in  a  decreased  effectiveness  but  still 
significantly  larger  than  its  OHD  counterpart. 
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CHAPTER  I 


STATEMENT  OF  THE  PROBLEM  AND  ASSOCIATED  LITERATURE 
1  .1  Introduction 

Magnetohydrodynamic  (MHD)  gas  flows  involve,  by  necessity, 
the  use  of  high  temperature  working  fluids  to  ensure  an  adequate  elec¬ 
trical  conductivity  even  with  the  addition  of  readily  ionizable  seed 
materials  such  as  the  alkali  metals.  In  particular,  for  MHD  generators, 
temperatures  of  approximately  2500°C  are  encountered  in  combustion  or 
open-cycle  equilibrium  generators  whereas  somewhat  lower  fluid 
temperatures ,  approximately  1800°C,  are  required  for  non-equilibrium, 
closed  cycle  generators.  Since  one  of  the  major  obstacles  to  achiev¬ 
ing  acceptable  generator  duty  cycles  between  maintenance  shutdown  is 

4. 

the  problem  [1]  of  materials  which  can  withstand  high  generator-duct 
temperatures, it  would  seem  meaningful  to  investigate  the  possibility 
of  protecting  the  MHD  duct  walls  from  these  high  temperatures  by 
cooling  techniques  which  have  proven  to  be  extremely  effective  in 
ordinary  hydrodynamics  (OHD)  flow  applications. 

Cooling  of  MHD  channel  walls  is  usually  accomplished 
externally  by  means  of  circulating  water  in  contact  with  the  exterior 
surface  of  the  walls,  by  water-jacket  designs  or  possibly,  when 
considering  the  insulator  walls,  by  a  specially  constructed  composite 
wall  designed  to  have  a  high  overall  thermal  conductivity  to  permit 
a  large  heat  transfer  rate  while  at  the  same  time  remaining  a  good 

^Numbers  in  square  brackets  indicate  references  listed. 
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electrical  insulator.  In  OHD  flows  it  is  known  that  cooling  methods 
involving  mass  transfer  at  the  surface  in  contact  with  the  hot  fluid 
stream  are  much  more  effective  than  cooling  the  wall  on  the  inside. 

It  was  therefore  suggested  to  investigate  the  possibility  of  transpi¬ 
ration  cooling  applied  to  the  insulator  or  B-walls  of  an  MHD  generator 
and  some  results  representing  first  order  effects  obtained  from  simple 
integral  methods  with  permeable-wall  boundary  conditions  are  presented 
as  a  separate  study  found  in  Appendix  A.  Objections  to  this  effort 
were  later  raised  when  considering  the  practical  applications  of 
constructing  such  a  porous  MHD  duct  and  a  film-cooling  technique  was 
instead  proposed  for  this  investigation. 

The  need  to  protect  surfaces  which  are  subjected  to  high 
temperature  environments  has  existed  for  other  than  MHD  flows  for  a 
long  time.  Generally,  the  high  temperature  environment  is  a  gaseous 
one  such  as  those  experienced  in  the  ionized  air  stream  surrounding 
re-entry  vehicles,  rockets,  plasma  jets,  and  high  performance  gas 
turbi nes . 

Sophisticated  methods  of  cooling  have  been  employed  over  the 
last  30  years  which  include  the  introduction  of  a  secondary  fluid 
into  the  boundary  layer  of  the  surface  which  must  be  protected. 
Ablation  is  another  method,  which  is  well  known  from  its  application 
to  re-entry  vehicles  in  NASA's  space  program.1  In  such  methods,  an 
added  coating  to  the  surface  known  as  a  "heat  shield"  decomposes  and, 
through  subl imination  and  other  highly  endothermic  processes,  a  large 
quantity  of  heat  is  absorbed  by  the  heat-shield  material.  Transpirati 


+Brewer,  R.  A.  and  Brant,  D.  N.,  "Thermal  Protection  System  for  the 
Galileo  Mission  Atmospheric  Entry  Probe"  AIAA  Paper  No.  80-0358,  1980 
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cooling  (ie.  cooling  Dy  mass  transfer  or  sweat)  is  a  second  method 
which  can  be  applied  to  a  surface  made  of  porous  materials.  The  cooler 
secondary  fluid  "diffuses"  through  such  a  permeable  surface  before 
entry  into  the  boundary  layer. 

Ablation  and  transpiration  cooling  are  primarily  effective 
in  protecting  the  region  where  the  materials  ablates  and  the  secondary 
fluid  enters  the  boundary  layer.  In  addition  to  heat  absorption  by 
this  material  or  fluid,  the  boundary  layer  is  effectively  thickened  by 
the  mass  addition  which  further  reduces  the  heat  transfer  rate  due  to 
smaller  surface  gradients.  The  short-comings  of  these  two  methods, 
however,  are  that  ablation  cooling  is  limited  to  systems  characterized 
by  high  heat  fluxes  of  short  duration  and  transpiration  cooling  re¬ 
quires  porous  materials  which  do  not  possess  the  necessary  structural 
strength,  which  in  most  cases  is  a  necessity  in  applications  where  such 
high  temperatures  are  encountered.  The  pores,  furthermore,  have  a 
tendency  of  becoming  clogged  resulting  in  relatively  short  duty  cycles. 

A  third  effective  method  of  cooling  involving  secondary  fuild 
injection  into  the  downstream  boundary  layer  is  termed  'film-cooling'. 
The  primary  objective  of  film-cooling  is  to  keep  the  surface  down¬ 
stream  of  the  injection  slot  cool.  Figure  1.1  represents  schematically 
a  typical  film-cooling  slot  which  would  be  installed  at  discrete 
locations  along  a  surface  to  protect  it  not  only  in  the  immediate 
region  of  the  fluid  injection,  but  also  for  a  considerable  distance 
downstream.  The  flow  far  downstream  of  the  coolant  injection  slot 
is,  in  fact,  similar  to  that  downstream  of  a  sweat  or  ablation- 

Furthermore,  film-cooling  could  be  used  to  shield 


cooled  surface. 
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a  surface  from  radiative  heat  transfer  in  addition  to  reducing  the 
convection  heat  transfer  rate  from  the  hot  main  stream  to  the  exposed 
wall.  Radiative  shielding  could  be  effected  with  gas  particle  suspen¬ 
sions  or  by  use  of  a  liquid  coolant.  The  present  work  does  not 
consider  such  processes,  however,  and  is  restricted  only  to  the 
influence  of  MHD  body  forces  and  plasma  non-equilibrium  on  the  film¬ 
cooling  effectiveness  which  will  be  defined  next.  In  addition,  only 
two-dimensional  compressible  flows,  including  axi symmetrical  geometries, 
will  be  considered  here. 


Fig.  1.1  Theoretical  Model  of  Film-Cooling  Slot 
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1.1.1  Adiabatic-Wall  Film-Cooling  Effectiveness 

For  the  thermal  design  of  a  surface  in  contact  with  a  high 
temperature  fluid,  the  engineer  requires  information  on  the  heat-trans¬ 
fer  rate  between  the  surface  and  fluid  and  the  variation  of  the  surface 
temperature  along  the  direction  of  flow  when  heat  transfer  at  the 
fluid/surface  interface  is  prevented.  This  would  represent  the  maxi¬ 
mum  temperature  to  which  the  surface  would  be  subjected  without  inter¬ 
nal  wall  cooling  and  is  known  as  the  adiabatic  wall  temperature,  1^. 

The  general  problem  in  film-cooling  theory  is  to  predict  the  relation¬ 
ship  between  the  wall  temperature  distribution  and  heat  transfer, 
given  a  certain  mainstream  flow,  secondary  injection  flow  and  specific  slot 
geometry.  The  design  problem  posed  may  also  be  to  predict  the  secon¬ 
dary  flow  needed  to  maintain  the  wall  surface  temperature  below  a 
certain  permissible  value  for  a  given  mainstream  flow,  slot  geometry, 
and  allowable  wall  heat  transfer.  The  heat  transfer  is  characteri zed 
by  the  Stanton  number  ST,  and  the  adiabatic  wall  temperature  by  a 
film-cooling  effectiveness,  n,  and  given  both  these  values,  the 
engineer  is  told  all  that  he  needs  to  know  for  further  design. 

The  Stanton  number  ST  in  the  presence  of  film-cooling  is 
defined  in  the  conventional  way,  Eq .  (1.1) 


or  since 


ST  =  q  /[(h  -h  )prUr] 
V  LV  r  w'  G  GJ 
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(i.i) 


(1.2) 
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and 


(1.3) 


ST  =  V[Cn(VU  +  7  Ur/] 


p  G  w' 
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(1.4) 


Equation  (1.4)  is  written  in  terms  of  quantities  which  the 

designer  can  easily  evaluate  together  with  what  he  desires  to  know, 

which  is  the  heat  flux,  q  ,  to  the  wall. 

'w 

The  adiabatic  wall  temperature,  T  ,  is  a  function  of  the 

aw 

primary  and  secondary  flow  fields,  slot  geometry,  and  temperatures  of 
the  two  fluid  streams.  By  defining  a  dimensionless  adiabatic  wall 
temperature,  q,  called  the  film-cooling  effectiveness,  this  temperature 
dependence  can  be  eliminated.  For  constant  property  flows  the  film¬ 
cooling  effectiveness  is  defined  by 


n 


aw  TG 

Tr  -  TG 

V-/  ^ 


(1.5) 


At  the  injection  slot,  T  =  Tr,  and  n  =  1 .  Generally 

aw  L 

T  <  Tr  and  Tr  <  T  so  that  n  varies  from  unity  at  x  =  0  to  0  at  x  =  °°, 
aw  u  b  g 

(see  Fig.  1.1)  where,  because  of  continual  mixing  of  the  two  gas  streams, 
the  adiabatic  wall  temperature  approaches  that  of  the  free  stream. 

For  variable  property  flows,  the  wall  temperature  in  the  slot 
would  be  expected  to  be  the  recovery  temperature  of  the  injected  gas 
whereas  far  downstream,  the  adiabatic  wall  temperature  would  be  the 
recovery  temperature,  T  ,  of  the  mainstream  gas.  As  a  consequence,  an 
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appropriate  expression  often  used  for  compressible  film-cooling  is 


(1.6) 


which  is  seen  to  be  equivalent  to  Eq.  (1.5) when  compressibility  effects 
can  be  neglected. 

For  the  large  temperature  differences  encountered  in  high 
temperature  MHD  channel  flows,  an  effectiveness  similar  to  Eq.  (1.6) 
will  be  employed  except  that  the  correspondi ng  enthalpies  will  be 
substituted  for  temperatures .  Hence  in  the  present  work,  the  film¬ 
cooling  effectiveness  is  defined  as 

h  -  h 


(1.7) 


which  can  be  rewritten  in  term  of  temperatures  as 

r^C2 

Taw  -  <TC  +  2  CT> 

n  = - F - j-  0-8) 

<VV  +  zh  (UG  "UC  > 

PG 

from  which  the  designer  can  determine  T  in  terms  of  known  quantities. 

Most  investigations  into  film-cooling  have  dealt  with  the  dete 
mination  of  heat  transfer  coefficient  and  adiabatic  wall  temperatures  as 
separate  studies  with  a  major  emphasis  on  predicting  and/or  measuring 
the  film-cooling  effectiveness,  q.  A  common  finding  is  that  the  heat 
transfer  coefficient  is  close  to  that  without  injection  which  implies 
that  it  is  primarily  determined  by  the  mainstream  flow.  The  adiabatic 
wall  temperature,  however,  varies  significantly  from  that  without 
secondary  gas  injection  which  makes  it  more  difficult,  and  perhaps 
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more  important,  to  predict.  In  the  present  work,  only  the  effects  of 
non-equilibrium  flow  and  MHD  body  forces  (both  on  average  velocities 
and  turbulent  fluctuations)  on  the  flow  and  temperature  fields  and 
adiabatic-wall  film-cooling  effectiveness,  will  be  investigated. 

1.1.2  Film-Cooling  Slot  Geometry  and  Theoretical  Flow  Model 

The  film-cooling  effectiveness  and  wall  heat  flux  are  depen¬ 
dent  on  a  large  number  of  variables  such  as  slot  geometry,  properties 
of  both  primary  and  secondary  flows,  initial  velocity  and  temperature 
(density)  ratios  of  the  two  gas  streams  (often  correlated  in  terms  of 

a  blowing  parameter  Q  =  p  Ur/prUr),  intensity  of  freestream  turbulence, 

CL  b  b 

etc . 

Experimental  investi gations.  requi re  diligent  efforts  to  ob¬ 
tain  many  data  records  and  lengthy  modifications  to  change  the  experi¬ 
mental  apparatus.  Often  the  results  from  experimental  measurements 
proved  the  effort  not  to  be  worth  while.  If,  on  the  other  hand,  a 
reliable  theoretical  prediction  method  can  be  devised  which  accurately 
models  the  film-cooling  process,  the  theoretician  needs  only  change  a 
few  statements  in  the  computer  program  and  the  results  are  relatively 
quickly  obtained. 

One  such  prediction  method  is  that  devised  by  Patankar  and 
Spalding  [2]  and  has  been  successfully  applied  to  OHD  film-cooling 
predictions  in  Ref.  [3].  Unlike  many  other  methods  of  predicting 
boundary  layer  flows,  that  of  Ref.  [3]  can  handle  velocity  profiles 
exhibiting  maxima.  It  has  therefore  also  been  adopted  in  the  present 
work  and  was  modified  in  a  general  way  to  allow  analysis  of  MHD  flows 


as  is  detailed  in  Chapter  IV.  The  theoretical  model  of  the  MHD  film¬ 
cooling  slot  is  shown  in  Fig.  1.1.  The  boundary  layer  flow  under 
these  conditions  will  be  turbulent  and  thermal  non-equilibrium  and 
finite  rate  ionization  will  be  considered. 

The  method  of  solution  is  of  the  finite-difference  variety 
and  one  main  feature  is  that  the  computational  grid  of  nodes  fits  the 
boundary  layer  exactly;  its  width  grows  or  contracts  at  the  same  rate 
as  the  boundary  layer  itself.  This  enables  the  same  number  of  nodes 
in  the  y  (or  w)  direction  to  accurately  calculate  dependent  variable 
gradients  as  the  marching  procedure  advances  downstream.  Furthermore, 
in  the  near  wall  regions  of  the  boundary  layer  where  gradients  change 
rapidly,  the  Couette  flow  analysis  of  the  present  method  precludes  the 
requirement  of  a  very  large  number  of  nodal  points  which  is  character¬ 
istic  of  other  methods  which  employ  constant  grid  spacings  across  the 
entire  width  of  the  boundary  layer.  For  such  methods  the  number  of 
nodal  points  is  large  by  necessity  and  must  increase  as  downstream 
solutions  develop  in  order  that  gradients  near  the  wall  may  be  calcu¬ 
lated  sufficiently  accurate.  In  the  present  method  the  length  of  com¬ 
puting  time  is  significantly  reduced  by  virtue  of  a  smaller  number  of 
grid  points  which  typically  may  only  be  of  the  order  of  50  while  other 
methods  require  of  the  order  of  1000  for  similar  accuracy. 
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1.2  Review  of  Associated  Literature 

1.2.1  Non-Equilibrium  MHD  Channel-Flow  Research  and  Developments 

The  effects  of  thermal  non-equilibrium  and  finite  rate  recom¬ 
bination  were  studied  in  the  insulator  boundary  layer  of  a  potassium- 
seeded  nitrogen  core-flow  MHD  Hall-neutralized  Faraday  accelerator  by 
D.  W.  Cott  [4,5].  The  non-similar,  compressible  boundary  layer  was 
assumed  to  be  laminar  and  two  dimensional.  The  general  multi-fluid 
conservation  equations  were  derived  from  Boltzmann's  equation  and  a 
two-fluid  model  was  adopted  consisting  of  the  over-all  gas  and  the 
electron  species.  A  collisionless  sheath  was  assumed  and  matched 
with  the  boundary  layer  equations  which  were  then  solved  numerically 
using  the  method  of  Ref.  [6].  For  the  laminar  flow  considered  little 
accuracy  was  lost  by  solving  the  linearized  equations  all  the  way  to 
the  wall  without  a  Couette-flow  model  for  the  overall  gas.  For  the 
electron  gas  boundary  layer,  however,  the  Couette-flow  concept  offered 
a  convenient  method  of  coupling  the  sheath  mixed  boundary  conditions 
to  the  linearized  electron  diffusion  and  energy  equations.  An  itera¬ 
tion  procedure  was  devised  which  matched  the  sheath  transition  region, 
treated  as  the  Couette-flow  region,  with  the  finite-difference  solution 
of  the  outer  edge  boundary  layer  flow.  It  was  concluded  that  thermal 
non-equilibrium  can  lead  to  significant  B-wall  shorting  in  long 
channels,  and  that  Hall  effects  should  not  be  neglected  when  large 
B-fields  are  encountered  and  that  operation  is  not  noticeably  affected 
by  ionizational  non-equilibrium  or  the  physics  of  the  electrostatic 


sheath . 
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High  and  Felderman  [7],  considered  the  compressible, 
turbulent,  boundary  layer  flow  over  the  B-wall  of  an  MHD  accelerator 
and  solved  the  equations  using  the  same  implicit  finite-difference 
method  of  Patankar  and  Spalding  [6].  The  multi-fluid  conservation 
equations  of  D.  Cott  [4,5],  were  formulated  for  two-dimensional 
unsteady  turbulent  flow  over  a  flat  plate.  The  flow  variables  were 
expressed  as  the  sum  of  a  mean  part  and  a  fluctuating  part  of  the 
form  Q  =  Q  +  Q'.  The  equations  were  rewri tten  i nto  divergence  form 
and  then  averaged  over  a  time  interval  that  is  large  compared  to  the 
fluctuations.  Following  classical  turbulent  boundary-1 ayer  theory 
only  pertubations  of  the  form  (pv)'Q'  were  retained. 

High  and  Felderman  modified  the  Couette  flow  curve  fits  of 
Ref.  [6],  given  for  the  global  momentum  and  energy  equations,  to 
include  the  MHD  terms  in  the  solutions.  For  the  electron  gas,  a  free 
molecular  sheath  was  also  assumed  to  exist  near  the  wall.  The  electron 
flux  from  the  continuum  analysis  is  matched  to  the  free-molecular  flux 
across  the  sheath  following  the  derivations  of  Sherman,  et  al .  [8]. 

The  resulting  expressions  used  were 
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The  continuum  expression  for  the  flux  on  the  left  hand  side 
contains  the  gradients  of  Cg  and  Tg  whereas  the  free  molecular  expres¬ 
sion  on  the  right  contains  the  values  of  Cg  and  T  .  Equating  the  con¬ 
tinuum  electron  energy  flux  at  the  edge  of  the  sheath  to  the  free- 
molecular  energy  flux  across  the  sheath  yields  a  second  relation 


I 


12 


k  3h 

C  <3T>w  +  I  <he+Ei> 

r 


3[C  (1+T  /!)■ 


By 


w 


ce  P  (f-) 

w  w 


3/2 


(-) 

7T 


1/2 


3/2  ™p  1/2  f 

ne  '8m. 
w  1 


Tim  1/2" 
2_£n4^8mT^ 


(1.10) 


where  again  a  relation  resulted  between  the  gradient  of  hg  and  the 

val ue  of  hg  at  the  wall, 
w 

For  the  Couette-flow  curve  fits  of  electron  energy  and  con¬ 
centration  the  wall  flux  could  be  expressed  in  terms  of  the  dependent 
variable  at  the  wall  by  using  the  boundary  conditions  of  Eqs.  (1.9) 
and  (1.10).  The  equations  were  then  solved  by  Runge-Kutta  methods  and 

curve  fit,  with  the  wall  values  as  parameters.  These  solutions, 

★ 

however,  assumed  a  constant  electron  current  jez  across  the  Couette 
layer  and  the  analysis  also  did  not  consider  any  MHD  damping  of  the 
turbulent  fluctuating  quantities. 

For  the  case  of  finite-rate  ionization,  a  =  a  (n  ,T  ),  no 
significant  amount  of  ionization  non-equilibrium  was  observed  in  the 
boundary  layer  except  near  the  wall  in  the  Couette-flow  region.  A 
small  amount  of  electron  thermal  non-equilibrium  was  observed  which, 
because  of  the  strong  dependence  of  the  electrical  conductivity  on 
electron  temperature,  caused  significant  differences  from  calculations 
using  thermal  equilibrium.  Inclusion  of  the  Hall  current  was  also 
shown  to  change  the  boundary  layer  characteristics  substantially. 

_  _  4* 

Sherman,  Yeh,  Reshotko,  and  McAssey  Jr.  [8,9]  studied  the 
compressible  boundary  layer  in  a  plasma  which  was  in  electron 

4- 

The  author  had  an  opportunity  to  discuss  various  aspects  of  this  study 
with  Dr.  E.  Reshotko  while  he  was  a  guest  speaker  at  the  University  of 
Alberta,  Mechanical  Engineering  Department. 
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thermal  non-equilibrium  and  in  which  ionization  and  recombination 
occurred  also  at  finite  rates.  The  surface  considered  was  partly 
an  insulator  and  partly  a  thermionical ly  emitting  electrode.  The 
flow  was  laminar  and  a  multi-fluid  MHD  model  was  assumed.  The 
governing  equations  were  solved  using  their  own  finite-difference 
scheme  and  results  obtained  showed  that  the  electron  temperature 
can  differ  significantly  from  the  heavy  particle  temperature.  The 
presence  of  finite  rates  influenced  the  electron  density  considerably . 
The  heat  transfer  due  to  electron  impact  and  electron  temperature 
gradient  was  considerably  less  than  the  heat  transfer  due  to  heavy 
particle  temperature  gradients  and  was  influenced  by  the  magnetic 
field  and  the  current. 

Other  related  studies  of  the  effects  of  thermal  non-equili¬ 
brium  and  finite-rate  ionization  are  contained  in  Refs.  [10],  [11],  and 
[12]. 

1.2.2  MHD  Turbulence  and  Turbulent  Flows 

A  relatively  recent  text  book  entitled  "MHD-Flows  and  Turbu¬ 
lence"  [13]  contains  two  excellent  sections  (Part  III  and  IV)  on  MHD 
Turbulence,  and  General  Problems  of  Turbulence.  Both  theoretical  and 
experimental  studies  are  presented. 

An  analytical  investigation  of  both  developing  and  fully 
developed  MHD  channel  flow  was  the  subject  of  a  PhD  thesis  by 
Scaglione  [14].  The  theoretical  derivation  of  Refs.  [15]  and  [16]  for 
a  two-fluid  model  was  utilized  to  describe  the  turbulent  MHD 
boundary  layer.  The  eddy  diffusion  of  momentum  in  the  near  wall 
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boundary  layer  was  based  on  a  modified  mixing-length  theory  according 
to  Refs.  [17],  [18]  and  [19].  In  the  outer  layer  the  eddy  viscosity 
was  based  on  a  modification  of  the  "intermittent  eddy-viscosity"  based 
on  studies  conducted  by  Klebanoff  [20].  Solutions  to  the  developing 
turbulent  MHD  boundary  layer  equations  were  arrived  at  by  a  linear¬ 
ization  technique  and  solved  with  an  implicit  finite-difference 

method  developed  by  Cebeci  and  Smith  [21,22].  Solutions  were  presented 

4  5 

for  channel  flow  Reynolds  numbers  ranging  from  10  to  4  x  10  and  for 
Hartmann  numbers  from  0  to  1200.  Results  showed  that  for  the  turbulent 
flow  case,  the  important  physical  parameters  such  as  displacement 
thickness,  momentum  thickness  and  skin-friction  coefficient  showed 
both  a  strong  Hartmann  and  Reynolds  number  dependence.  The  results 
also  exhibited  the  effect  of  detransition  from  turbulent  to  laminar 
flow  when  the  ratio  of  Hartmann  number  to  Reynolds  number  exceeded 
1/225.  Asymptotic  solutions  for  fully  developed  flow  showed  quite 
close  agreement  with  experimental  data  and  with  similar  theoretical 
analysis  of  Lykoudis  and  Brouilette  [17,18].  A  correlation  of  the 
skin-friction  coefficient  based  on  a  laminar  sub-layer  and  turbulent 
core  model  was  derived  with  Reynolds  and  Hartmann  numbers  as  parameters. 

This  agreed  to  within  10%  of  the  experimental  data  of  Ref.  [17].  The 
analysis  was  limited  by  the  assumptions  of  incompressible  flow, 
electrically  non-conducting  channel  walls,  small  magnetic  Reynolds 
number,  absence  of  an  externally  applied  electric  field,  a  magnetic 
field  fixed  relative  to  the  channel,  negligible  Hall  effects, 
electrically  neutral  fluid,  steady  state  flow,  and  two-dimensional 
constant-property  flow  in  which  only  electromagnetic  body  forces  are  present. 
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Experimental  measurements  of  turbulent  friction  factors  of 
electrically  conducting  liquids  in  pipes  and  the  influence  of  a  long¬ 
itudinally  applied  magnetic  field  were  reported  by  Levin  and  Chinenkov 
[23].  This  was  paralleled  with  a  theoretical  study  by  Levin  [24], 
which  adopted  a  similar  approach  taken  by  Lykoudis  and  Brouilette  [18] 
for  deriving  a  first  order  effect  of  the  MHD  body  forces  on  the  turbu¬ 
lent  velocity  fluctuations.  Experimental  data  was  obtained  for  a  wide 
range  of  Reynolds  and  Hartmann  numbers  and  results  showed  that  a 
substantial  increase  in  the  Hartmann  number  compared  to  those  reported 
in  the  literature  was  required  to  reduce  the  friction  factor  by 
approximately  an  order  of  magnitude. 

Other  relative  theoretical  studies  of  MHD  turbulent  flows 
and  MHD  turbulence  modeling  including  mixing-length  theories  are 
found  in  Ref  [25],  which  consider  compressible  turbulent  MHD  boundary 
layers.  The  work  in  Ref.  [26]  showed  that  the  power  output  of  a  Hall 
generator  depends  strongly  on  the  wall  temperature  under  laminar,  fully 
developed  flow  conditions,  but  that  this  dependence  is  greatly  reduced 
when  the  turbulent  boundary  layer  velocity  profile  was  assumed.  This 
would  imply  a  favourable  effect  of  film-cooling  on  MHD  generator  per¬ 
formance,  since  the  injection  of  a  secondary  coolant  fluid  results  in 
a  turbulent  boundary  layer  downstream  from  the  slot.  Moreau  [27], 
applies  spectrum  tensor  analysis  and  derives  spectrum  functions  and 
distribution  and  transfer  of  energy  in  wave-number  space  for  Joule 
dissipation  and  induced  electric  field.  The  study  is  restricted  to 
the  domain  of  small  magnetic  Reynolds  numbers,  Rem  «  1 ,  so  that  the 
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fluctuations  of  the  magnetic  field  are  negligible  compared  with  those  of 
the  velocity  since  |b|  ~  Rem  |u'|  /u.  A  system  of  equations  is  presented 
in  Refs.  [28,29]  for  calculating  the  average  and  pulsative  velocity  and 
temperature  characteristics  in  a  turbulent  jet  of  an  incompressible 
conducting  liquid  in  the  presence  of  a  longitudinal  or  transverse 
magnetic  field  (Re  <<1).  The  semi-empirical  Kolmogorov  relation  for 
the  energy  dissipation  and  the  Rotta  relation  for  the  exchange  term  were 
employed  as  well  as  employing  the  gradient  representation  of  the  turbu¬ 
lent  diffusion.  An  MHD  analogue  of  the  Prandtl  formula  for  the  friction 
stress  and  the  heat  flux  were  derived.  The  so-called  first  Prandtl 
formula,  for  example,  was  written  as 


T_ 

P 


-  r3/2C  2 

uv  =  e  ,  L  i 

u  k3 


9u  |  9u 

9y*  3y 


(1.11) 


where  is  a  function  of  the  local  interaction  parameter, 
a(o82/p)/(3u/9y) ,  with  a  =  1  -  6,  where  0  is  assumed  to  be  constant 
(0  <  3  < 1 )  and  must  be  determined  from  experiment. 

It  is  also  shown  that  for  some  limiting  value  of  the  local 
interaction  parameter  (Stuart  number)  there  is  complete  suppression  of 
the  turbulent  motion  which  in  turn  precludes  turbulent  friction.  With 
an  increase  in  Stuart  number,  the  MHD  Prandtl  number  increases  and  in 
the  limiting  case  approaches  unity. 

Further  experimental  investigations  into  the  effect  of 
electromagnetic  fields  on  turbulence  fluctuations  are  presented  in 
Refs.  [30,  [31],  and  [32]. 
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1.2.3  OHD  Film-Cooling  Research  and  Developments 

An  excellent  review  of  film-cooling  investigations  is  pre¬ 
sented  by  Goldstein  [33].  The  general  theory  and  method  of  analysis 
is  discussed  for  two-  and  three-dimensional  film-cooling  of  both 
incompressible  and  compressible  flows.  Also  included  is  a  summary  of 
run  conditions  and  geometry  details  for  40  experimental  investigations 
including  those  of  Wilson  [34,35].  The  influence  of  variable  free 
stream  velocity,  free  stream  turbulence,  film-cooling  slot  geometry 
and  foreign  gas  injection  on  the  film-cooling  effectiveness  and  heat 
transfer  are  discussed.  An  in-depth  review  of  three-dimensional  film¬ 
cooling  is  also  included. 

A  general  review  of  the  Patankar  and  Spalding  numerical 
solution  procedure  [2],  is  presented  by  Spalding  in  Ref.  [36].  Both 
the  general  behaviour  of  the  film-cooling  effectiveness,  n,  and 
Stanton  number,  ST,  are  analyzed  in  reference  to  variations  in  para¬ 
meters  such  as  slot  geometry,  secondary-to-primary  gas  velocity  ratios, 
free-stream  Mach  number,  coolant  composition,  etc.  The  general  method 
of  solution  evolved  from  earlier  film-cooling  investigations  at  Imperial 
College,  London,  such  as  those  conducted  by  Cole  and  Spalding  [3]. 

The  method  has  been  compared  favourabl  y  wi th  experimental  measurements 
involving  different  U^/Ug,  T^/T^,  pc/pG  ratl0S  as  we^  as  foreign  gas 
injection.  Allowance  for  free-stream  turbulence  is  incorporated  into 
the  numerical  finite-difference  procedure.  Further  refinements  were 
necessary,  however,  for  the  influence  of  slot  lip  thickness  for  which 
experimental  measurements  did  not  support  theoretical  predictions. 
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1.2.4  MHD  Film-Cooling  Research  and  Developments 

It  was  mentioned  earlier  that  one  potential  application  of 
film-cooling  techniques  is  for  plasma  jets  and  torches.  Shiver 
[37,38]  considered  film-cooling  a  plasma  torch  as  a  means  of  improving 
the  overall  energy  conversion  efficiency  of  the  torch.  Regenerative 
cooling  methods  had  been  successfully  used  to  cool  plasma  torches  and 
have  produced  higher  operating  efficiency,  defined  as  the  ratio  of  the 
energy  content  of  the  emerging  gas  stream  to  the  electrical  energy 
input.  Torches  with  transpiration-cooled  anodes  have  produced  effi¬ 
ciencies  of  80  to  90  percent  compared  to  50  to  70  percent  efficiencies 
for  the  conventional  water-cooled  torch.  The  work  by  Shiver  represented 
the  first  experimental  investigation  of  film-cooling  applied  to  plasma 
torches.  A  gas  stabilized  constricted-arc  type  torch  was  constructed 
with  interchangeable  water-cooled  copper  and  film-cooled  carbon  anodes 
for  the  use  with  argon.  A  shell  and  tube  calorimeter  was  constructed 
to  determine  the  energy  content  of  the  torch  effluent.  Tests  were  con¬ 
ducted  with  each  torch  in  which  arc  current  and  gas  flow  rate  were 
selected  as  the  operating  variables.  Two  test  runs  were  also  made 
using  the  film-cooled  torch,  during  which  the  arc  current  and  primary 
gas  flow  rate  were  maintained  constant,  while  the  secondary  (coolant) 
flow  was  varied. 

The  results  indicated  that,  for  approximatel y  the  same  total 
gas  flow  rate  and  arc  current,  the  torch  equipped  with  a  film-cooled 
carbon  anode  was  approximately  20  percent  more  efficient  than  a  fully 
water-cooled  torch.  The  film-cooled  torch  efficiency  ranged  from  80 
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to  90  percent,  while  the  water-cooled  efficiency  ranged  from  60  to  70 
percent  for  gas  flow  rates  between  0.20  Ibs/min.  to  0.57  Ibs/min.  An 
error  analysis  indicated  that  results  obtained  were  correct  to  within 

-  u. 

Theoretical  studies  of  either  transpiration  or  film-cooling 
applied  to  MHD  flows  include  those  of  D'Sa  [39],  who  considered  a 
hydro-magnetic  wall  jet.  The  effect  of  magnetic  forces  on  the  laminar 
wall  jet  of  an  electrically  conducting,  incompressible  fluid  of  constant 
properties  was  studied  theoretically.  The  imposed  magnetic  field  was 
assumed  to  be  uniform  and  normal  to  the  wall,  which  could  be  porous  or 
slotted  to  allow  suction,  blowing  or  tangential  injection.  The  analysis 
was  limited  to  laminar,  incompressible,  constant  property  flows.  The 
Levy-Lees  transformation  was  used  and  a  series  expansion  was  sought 
about  the  non-magnetic  solution  of  the  transformed  boundary  layer 
equations.  The  required  series  was  constructed  in  powers  of  the  so 
called  interaction  parameter,  N  ,  sometimes  referred  to  as  the  Stuart 

X 

number,  which  physically  represents  the  ratio  of  magnetic  to  inertial 
forces.  The  resulting  linear  equations  were  solved  using  a  modified 
Runge-Kutta  method  together  with  Newton-Raphson  iteration  for  three 
values  of  the  wall  blowing  (suction)  parameter,  representing  either  an 
impermeable  wall,  suction  or  injection. 

Luft  and  Rodkiewicz  [40]  applied  the  Karman- Pohl  hausen 
technique  to  obtain  solutions  of  the  MHD  boundary  layer  flow  over  an 
insulated  flat  plate  in  the  presence  of  a  uniform  transverse  magnetic 
field  and  uniformly  applied  wall  suction.  Incompressible,  constant 
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property  flow  was  also  assumed  and  a  further  transformation  of 
variables  following  the  method  of  Lew  (see  Ref. [9]  of  Appendix  A)  was 
performed  before  deriving  the  integral  equation  for  MHD  flate-plate 
suction  flow.  Asymptotic  MHD-suction  profiles  were  obtained  and  fully 
developed  boundary  layers  illustrated  for  combined  Hartmann  numbers  and 
suction  coefficients.  Developing  velocity  profiles  and  average  skin 
friction  coefficients  were  derived. 

The  similarities  observed  in  this  investigation  between  suc¬ 
tion  and  MHD  velocity  profiles  led  to  the  proposed  modification  of  the 
constant  A*  appearing  in  the  laminar  sublayer  portion  of  the  two-layer 
eddy  viscosity  model  of  turbulence  as  discussed  in  Section  3.2  of  this 
thesis.  The  complete  first  order  analysis  of  MHD  suction  boundary 
layer  flow  by  Luft  and  Rodkiewicz  is  presented  in  Appendix  A  and  is 
included  as  part  of  the  present  thesis  work. 

Sych  [41]  used  a  variational  method  to  investigate  the  steady- 
state  flow  of  a  conducting  fluid  in  a  straight  channel  of  rectangular 
cross  section  with  transverse  magnetic  field  when  there  is  forced  flow 
through  the  side  walls.  The  liquid  viscosity  was  taken  into  account 
as  well  as  the  conductivity  of  the  channel  walls  perpendicular  to  the 
external  magnetic  field.  Results  obtained  were  used  to  assess  the 
losses  in  an  MHD  generator. 

Jain  and  Mehta  [42]  presented  an  exact  solution  in  a  closed 
form  for  the  hydromagnetic  equations  for  an  incompressible,  viscous, 
and  electrically  conducting  fluid  flow  through  an  annulus  with  porous 
walls  in  the  presence  of  a  radial  transverse  magnetic  field.  Solutions 
were  reduced  to  a  simpler  form  by  introducing  the  normal  assumption  of 


I 


I 
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v/V|_|  <<  1  ie.  low  Magnetic  Reynolds  number,  and  the  pressure  distri¬ 
bution  was  discussed  for  this  case.  The  analysis  was  also  extended  to 
include  parallel  MHD  Couette  flow. 

The  present  work  investigates  the  adiabatic  film-cooling 
effectiveness  of  film-cooling  the  insulator  (B)  wall  of  both  an  MHD 
equilibrium  Faraday  generator  and  an  MHD  non-equilibrium  Faraday  gen¬ 
erator.  Turbulent  2-D  flow,  including  the  damping  of  turbulent  fluc¬ 
tuations  by  the  electromagnetic  fields  is  considered  allowing  for  plasma 
electron  thermal  non-equilibrium  and  finite-rate  ionization.  The 
implicit  finite-difference  method  of  Patankar  and  Spalding  [2]  is 
modified  in  a  general  way  to  allow  MHD  flow  calculations.  On  the  basis 
of  the  literature  review  presented  here,  such  a  film-cooling  investi¬ 
gation  has  not  been  undertaken  before.  Additional  research  papers 
consulted  in  preparation  of  this  thesis  are  found  in  Refs.  [43]  to 
[58]  inclusive. 


CHAPTER  II 


THEORETICAL  DEVELOPMENT 
2.1  Boundary  Layer  Equations 

The  development  of  the  equations  of  motion  for  a  two-dimen¬ 
sional,  turbulent,  non-equilibrium  plasma  flow  in  an  MHD  duct  follows 
closely  the  derivation  of  their  laminar  counterpart  given  in  detail  in 
Ref. [4].  The  classical  derivation  is  performed  where  the  variables  are 
expressed  as  the  sum  of  a  mean  part  and  a  fluctuating  part  after  which 
the  equations  are  averaged  over  a  time  interval  which  is  large  com¬ 
pared  to  the  period  of  the  fluctuations.  The  mixing-length  concept  is 
used  to  affect  closure  of  the  system  of  equations  and  only  the  boun¬ 
dary  layer  on  the  B-wall  (insulator  wall)  is  considered.  With  reference 
to  the  coordinate  system  shown  in  Fig  2.1,  the  following  equations  can 
be  shown  to  be  identical  with  those  given  in  Ref.  [7]. 


Fig.  2.1  Schematic  of  MHD-Channel  Co-ordinate  System 
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Equati ons  (2. 1)  to  (2 . 3) represent  the  conservation  equations  of 
the  heavy  species,  namely,  the  nitrogen  carrier  gas  and  the  potassium 
seed  neutrals  and  ions.  In  addition  to  the  conventional  boundary 
layer  terms  the  equations  include  terms  which  represent  the  MHD  body 
force,  Joule  dissipation,  transport  of  energy  by  ion  and  electron 
di ffusion, and  heat  conduction  due  to  the  presence  of  electron  temper¬ 
ature  gradients. 

Equations  (2.4)  and  (2.5)  describe  the  "electron  gas"  and 
included  are  terms  which  represent  finite- rate  ionization-recombination 
electron  diffusion,  energy  transport  due  to  electron  diffusion,  elec¬ 
tron  temperature  gradients,  electron  production,  Joule  dissipation  and 
collisional  energy  losses. 

When  considering  equilibrium  ionization,  Eq.  (2.4) is  not 
used  and  the  electron  mass  fraction,  C  , i s  computed  from  the  Saha 
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and  the  partition  function  can  be  approximated  by 


Zf$  =  9.115  x  10"1  +  7.33  x  10“5T  -  1 .49  x  10'8T2 

(T  <  6000°K) 


(2.10) 
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For  the  case  of  finite-rate  or  frozen  ionization  the  mass 
fraction  is  obtained  from  the  solution  of  Eq.  (2.5)  and  the  electron 
production  term  wg  is  approximated  by  an  empirical  expression  given 
in  Ref.  [5]. 
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2.2  Molecular  Transport  Properties 

The  molecular  transport  properties  required  in  the  governing 
equations  are  the  viscosity  of  the  gas,  the  electron  thermal  conduc¬ 
tivity  and  the  electrical  conductivity  of  the  gas.  Sutherland's  vis¬ 
cosity  law  is  used  to  prescribe  the  laminar  viscosity  of  the  carrier 
gas 

y  =  1.086  x  10~6  T3/2/(T  +  100°K)  N-sec/m2 

Other  properties  have  been  taken  from  Ref.  [5].  The 
electron  thermal  conductivity  is  obtained  from 
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where  and  Tq  are  defined  in  Refs.  [4,59].  The  electrical  conductivity 
is  given  by  the  expression 
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with  the  el ectron-el ectron  interaction  correction  given  by 
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where  the  expressions  for  A^  ,  A^  ,  A  ,  Q  •  are  defined  in  Ref. [60] 
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The  equation  of  state  as  used  in  Ref. [5]  is 
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where  T  is  in  °K  and  p  is  in  atmospheres. 

The  collisional  energy  term  in  Eq-  (2.5) (last  term)  is  derived 

6 

from  the  expression  3kC  p  S  .1  v*  (T  -T  )  as  is  shown  in  detail  in 

e  s^e  ms  es  s  e 

Ref.  [4]  where  is  the  fractional  energy  loss  during  an  inelastic 

collision  of  an  electron  with  specie  "s"  (6  ,  . .  =1)  and  these 

elastic 

inelastic  losses  are  lumped  into  one  term,  6  defined  by  Ref.  [61]. 
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Experiments  performed  in  Ref. L 62] i ndi cate  that  for  the  nitro¬ 
gen  background  gas  (6^  =  18)  considered  here,  6  ^  remains  more  or 

2_4 

less  constant  at  7  X  10  for  temperatures  in  the  range  of  2800  K  to 


5000 


-3 


°K.  In  contrast,  for  air,  Ref. [7]  gives  6  ^  =  2.93  X  10  which 


corresponds  to  a  6  .  of  about  150.  Needless  to  say,  the  properties  of 

air 

actual  combustion  gases  in  real  MHD  generator  installations  must  be 
accurately  established  before  a  meaningful  comparison  between  the  theory 
and  experiment  can  be  made. 

The  Prandtl  and  Schmidt  numbers  were  assumed  to  be  equal 
with  the  laminar  values  of  0.7  and  their  turbulent  counterpart  assigned 
the  values  of  0.9. 


2.3  Boundary  Conditions 

The  dependent  variables  in  Eqs.  (2. 1)  to  (2.6)  must  be  specified 
at  the  channel  insulator  wall,  in  the  free  stream,  and  at  the  film¬ 
cooling  slot  exit  to  begin  the  integration. 


Boundary  values  at  the  wall  for  velocity  and  enthalpy  are: 


.  OH 

u=v=o  and 


w  =  0.  For  the  electron  mass,  C  ,  an  electron  sheath 

’  e,w 


analysis  is  performed  following  the  procedure  of  Sherman,  et  al . 

[8,9].  A  collisionless  free  molecular  sheath  is  assumed  to  exist 
at  the  wall  and  the  electron  flux  across  the  sheath  is  matched  with 
the  continuum  electron  flux  at  the  edge  of  the  sheath  resulting  in  [7] 
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A  similar  analysis  for  the  electron  energy  flux  results  in 
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Equations  (2.19)  and  (2.20)  serve  as  mixed  boundary  conditions  for  the 

electron  mass  fraction  C  and  electron  enthalpy  h  ,  respectively,  for 

e  e 

the  case  of  frozen  and  finite-rate  ionization-recombination.  When 

equilibrium  ionization-recombination  is  considered,  Eq.  (2.20)  only 

is  used  as  boundary  condition  for  h  from  which  C  can  be  evaluated 

e  e 

using  the  Saha  equilibrium  relationship. 

The  electric  and  magnetic  field  parameters  E^  and  B  are 
specified  independently  and  are  assumed  constant  across  the  generator 
duct.  Actually  Ez  is  calculated  from  a  given  load  parameter  (often 
referred  to  as  either  a  generator  or  a  pump  coefficient)  defined  by 


K  _  +  z  _  applied  electric  field 
1  —  BIL  induced  electric  field 


(2.21) 
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where  the  sign  of  is  positive  when  Ez  is  in  the  same  direction  as 
the  induced  field  (accelerator  mode)  and  negative  when  Ez  is  directed 
opposite  to  BUg  as  in  the  case  of  a  generator.  According  to  a  study 
by  Swift-Hook  and  Wright  [63],  an  optimum  generator  design  is 
achieved  when  | Kj  -  .75  and  this  is  the  value  prescribed  at  the 
film-cooling  slot  exit.  With  an  appropriate  variation  in  the  external 
load  resistance,  Ez  can  be  kept  constant  by  allowing  to  vary  as 
with  downstream  distance, x. 

The  effects  of  Hall  currents  can  be  assessed  by  setting 
E  =0  and  by  applying  an  electric  field  E  such  that  the  free  stream 

X  X 

only  is  Hall  neutralized  (je  =0)  with  E  calculated  from  Eq.  (2.6) 

X  5  vjj  X 

evaluated  in  the  free  stream.  When  Hall  currents  are  precluded  or 
neglected  then  je  =0  everywhere. 

X 

The  free  stream  boundary  conditions  are  specified  from  the 
solution  of  the  one-dimensional  core-flow  equations  which  may  be 
obtained  from  Eqs  .  (2 . 1)  to  (2 . 6)  by  dropping  all  transverse  gradient 
terms.  Furthermore,  previous  investigation  of  non-equilibrium  plasma 

channel  flows  [5,7],  indicate  that  the  core  qas  is  essentially  in 
ionization-recombination  equilibrium  so  that  the  electron  number 
density  or  mass  fraction  can  be  calculated  from  the  Saha  equilibrium 
relations,  Eqs. (2 . 7)  to  (2 . 10),  based  on  the  free  stream  electron  temper¬ 
ature, T  . 
eG 

The  resulting  equations  are  then 


prlLA  =  m  =  constant  mass  flow 
b  b 


(2.22) 
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dU, 


PrU 


G  G  dx 


=  '  (l^+  Jex,GB) 


(2.23) 


dHp 

PGUG  dx-  =  jez,GEz  +  ^ex ,GEx 


(2.24) 


prUP  ~  (C*h  ) 
G  G  dx  e  e 


P  ke2Pf2pCer 

Jez.GEz  +  Jex.GEx  +  f-^^  WT-V  <2'25) 


■x,G  x  2  3  o0 


jez,G  =  °e,GUGB(1+KL) 


(2.26) 


pe  =  Zpgrtg 


(2.27) 


Since  there  are  only  six  equations  with  nine  unknowns, 

namely  pg,  Pg,  Ug,  Hg,  jez,g,  hg  g,  A,  B,  and  Ez  (or  )  three  of 

these  must  be  specified  independently.  Cott  [4,5]  specifies  Tr, 

b 

je  r,  and  B  whereas  High  and  Felderman  [7]  specify  hp,  E  ,  and  B. 

Z  ,  b  G  Z 

A  constant  area  duct  is  considered  here  so  that  A,  B,  and  E 
(or  K^)  have  been  assigned  independently. 

Numerous  solutions  for  one-dimensional  constant  area  non¬ 
equilibrium  generators  have  been  presented  by  Heighway  and  Nichols, 

[10]  for  both  supersonic  and  subsonic  velocities  and  a  typical 
solution  providing  the  outer  boundary  conditions  is  depicted  in 

Fig  2.2. 
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Fig.  2.2  Typical  MHD1-D  Core  Flow  Solutions,  Equilibrium  Case 
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2.4  Calculation  of  Axial  Flow  Variables 

Downstream  pressure  is  calculated  in  accordance  with  Patankar 
and  Spalding's  procedure  for  handling  confined  flows  [2]  using  the 
pi ug-f 1 ow  formul a  ^ 


F' 


dx 


1 


2u_  dm  mu  dT 
A  dx  at  dx 

-  mu/Ap 


(2.28) 


where  T  is  an  average  temperature,  u  is  the  mass-average  velocity 
defined  by 


u 


pu2dy/  _ 

duct  duct 


pudy 


and  the  retarding  force  per  unit  duct  length  F'/A  is  given  by 


+  je^B 


Once  the  pressure  gradient  has  been  obtained  from  Eq - (2.28) 
the  downstream  local  enthalpy  in  the  core  flow  may  then  be  evaluated 
from 


dh 


G 


dx 


_L  +  J6z,G 
Pq  Pq 


(Kl+1)B 


(2.29) 


^  A  misprint  in  Ref.  [2]  is.noted  here  in  the  denominator  which  is 
incorrectly  written  as  1 -mu/Ap 


Previous  investigations  by  Cott  [4,5]  and  High  and 
Felderman  [7],  showed  that  ionizational  equilibrium  exists  in  the 
core  flow  so  that  the  downstream  electron  enthalpy  can  be  determined 
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from 

dh 


dx 


K.  Bje  r  eB(K,  +1 )  ~ 

1  z’b  +  L  o+3 

prC  m  PG  2 

be-  e 


ke2 

m  3a 
e  o 


6effpGC( 


(VV 


(2.30) 


and  the  downstream  electron  mass  fraction  can  then  be  computed  from 
the  Saha  relation  using  the  electron  enthalpy  computed  from  Eq.  (2.30). 


2.5  Initial  Conditions 

The  velocity  profile  shown  in  Fig.  2.3,  at  the  slot  exit  is 
based  on  a  power-law  variation  of  u  with  y  and  is  one  of  several 
profiles  considered  by  Cole  and  Spalding  [3]  in  their  work  to 
assess,  along  with  other  parameters,  the  influence  of  different 
initial  profiles  on  the  theoretical  prediction  of  the  film-cooling 
effectiveness.  It  is  essentially  a  power-law  profile  prescribed  as 
fol 1 ows : 

For  6w  u/Uc  =  (y/6w)  W  (2.31 ) 


For  6^ 


For  6^ 


A 


(u-UL)/(Uc-UL) 

=  (yc-y)/\ 

(2.32) 

(u-UL)/(UG-UL) 

r  i nG 

=  (y_>c"V/<sG 

(2.33) 
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To  define  the  initial  velocity  profile,  therefore,  the  quan¬ 
tities  Ur,  U,  ,  Up,  6.  ,6p,  6  ,  n  ,  n,  ,  and  n.  must  be  specified.  The 
tLGLGwwL  G 

initial  thickness  of  the  boundary  layer,  y  ,  identifies  the  grid  line 
go  =  1  ,  which  defines  the  outer  edge  of  the  layer  as  it  develops  down¬ 
stream.  The  inner  edge  of  the  layer  at  y  =  0,  coincides  with  the  film- 
cooled  wall  where  the  grid  line  has  the  value  oo  =  0. 

Initial  gas  and  electron  temperature  profiles  were  assumed 
equal  and  uniform  across  the  slot  and  in  the  free  stream  and  varying 
linearly  across  the  slot  lip.  Hence  for 

0  1  y  <  :  Tc  =  constant  (2.34) 

ycU<  (yc+tL)  :  (T-Tc)/(Tg-Tc)  =  (y-yc)/tL  (2.35) 

y  >_  (y£+t|_)  :  Tg  =  constant  (2.36) 


Finally,  the  electron  mass  fraction,  C  ,  was  calculated 
assuming  ionizational  equilibrium  from  Eq.  (2.7).  Thermal  and  ion- 
izational  equilibrium  are  assumed  in  defining  initial  boundary  layer 
profiles  for  both  equilibrium  and  non-equilibrium  MHD  flows  considered 
in  this  work. 

Because  of  numerical  computation  problems  experienced  by 
assuming  a  step  function  for  the  magnetic  field  parameter  B  (from  0  to 
B  =  constant  at  x/yc  =  0),  the  initial  value  for  B  was  applied  as  a 


linear  function  of  x  until  a  finite  value  of  x/y^  was  reached,  after 
which  B  remained  constant  with  x.  The  effect  of  choosing  different 
values  of  x/y^  at  which  B  becomes  a  constant  is  assessed  and  dis¬ 
cussed  in  Section  5.3.4.  Typical  values  of  x/y^  assigned  were  20,  10, 
and  1,  where  the  smaller  value  of  x/y^  approaches  closer  to  a  step 
function  for  B(x). 


Fig.  2.3  Initial  Velocity  Flow  Model  Assumed  in  Present  Work 


CHAPTER  III 


SOLUTIONS  TO  THE  GOVERNING  EQUATIONS 


3.1  Numerical  Solution  Procedure 

The  numerical  method  adopted  here  to  solve  the  MHD  boundary 
layer  equations  is  that  of  Patankar  and  Spalding  first  described  in 
detail  in  Refs.  [2]  and  [6].  It's  particular  application  to  the  OHD 
film-cooling  process  is  described  in  generality  in  Ref.  [64]  and  a 
specific  application  to  evaluate  the  effects  of  such  parameters  as  slot 
lip  thickness,  free-stream  turbulence,  initial  velocity  profiles,  free- 
stream  to  slot  velocity  ratios  etc.  on  theoretical  film-cooling  pre¬ 
dictions  is  demonstrated  in  Ref.  [3]. 

The  Patankar-Spalding  formulation  transforms  the  governing 
equations  from  the  Cartesian  to  the  V on  Mises  coordinate  system  with 
the  change  of  variables  defined  by: 


£  =  x,  \p  =  /pudy 

such  that 

JL  =  i_vJL_3_=nlljL  (3-1  ) 

9x  £  ~  u  9 y  ’  9y  p  9^ 


The  transformed  momentum  equation  (2.2),  for  example,  is 
then  readily  developed  by  substitution  of  (3.1)  to  obtain 


pu 


v_  _9u 
u  9y 


3u  _ 
+  pv  —  = 
9y 


9p  v  9p 
9C  "  u  9y 


je  B 
z  y 


+ 


Vfpu 


9u 

9^ 


(3.2) 


36 


37 


Replacing  E,  by  x  after  having  performed  the  differentiation  and 
rearranging,  Eq .  (3.2)  becomes 


3u 

3x 


_a_ 

dip 


Meff  Pu 


3u 

3T> 


.  J_  (|E+  je  B  ) 
pu  v3x  J  z  y' 


(3.3) 


The  x  -  i|i  grid  is  next  confined  to  the  boundary  layer  region 
only,  by  a  non-dimensional izinq  of  Eq.(3.3)  with  respect  to  ip  by 
defining  a  dimensionless  stream  function 


ip  -  ipj 


(3.4) 


where  ^  and  are  the  values  of  the  stream  functions  at  the  inner  and 
outer  boundaries  respectively.  The  effect  is  a  grid  which  grows  along 
with  the  boundary  layer  in  a  continuous  manner  and  a  cross-stream 
variable  which  varies  always  between  zero  and  unity. 

The  momentum  equation  (3.3)  in  terms  of  x  and  co  then  becomes 


3u 

f  ^  | 

mi 

'  +u)(m£  '  -mj 1 

')  ' 

3u 

3 

[  yeff  pu 

> 

3u 

1 

3x 

l 

^ j 

> 

3u> 

3oj 

(>pE-<pj)2 

3oo 

j 

PU 

3x 


’z  y‘ 


(3.5) 

For  the  sake  of  completeness,  the  remaining  boundary  layer 
equations  (i.e.  Eqs.  (2.3)  through  (2.5))in  terms  of  x  and  co  are  written 
as  follows: 


Overall  Energy 


3H  , 

'  mj 1  +co(m^ 1  -mj 1  ^ 

3H  _  3 

pu  peff ,H  3H 

3?  + 

^E  -  ^ 

3a)  3u) 

/  .  ,  \2  3oo 

'  °eff ,H  ' 

38 


9u) 


puy 


eff ,  H 


1- 


1 


eff,H 


a  ii2  Pu  hn 
~  (— )  +  e 

9cu  2  +  7  7 


y 


eff  ,C 


a 


eff  ,C 


eff  ,H 
Wf,H 


9C 
_ ( 

9co 


+  ^ 


pu  h 


A  (!e. 

V  T  > 


s(vy)2  3^T 


je  E  +  je  E 
XX  0  z  z 

~pu 


J 


(3.6) 


where 


y 


eff  ,C 
7eff  ,C 


e_  _  y 
S 


1  + 


T 


(3.7) 


Meff  ,H  =  jj_  + 

aftrr  .  P  P  + 
eff,H  r  t 


(3.8) 


Electron  Conservation 


5 

9x 


mj  '+co(m^  '  -mj '  j 


^E  ■  ^ j 


!^e=  j.  f 

9oo  9oo 


pu 


y 


eff ,C  9C  ^ 

e  e 


(^E-^j)2  aeff ,ce 


9co 


w 

+  —  +  — 

pu  9co 


pu 


.(^E"^j ) 


.  T  1 
2  ve  S  9w  1  T  j 


(3.9) 


Electron  Energy 


9C  h 
e  e 
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9C  h 
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(^E-^j) 


2  9co 
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e  (^E-^j) 


Vi  +  EzJez  *  Exjex  +  eCe 


PU 
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-  E 
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-  2  2 
^  k  0  PC  (5  rr 
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(3.10) 
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where 


y 


eff,e  e  .  e 

j  ,,  Cn  C  P+ 

eff.e  P  e  t 


(3.11) 


These  equations  share  the  common  form 


|i  +  (a+bu)  |i  =  (c  +  d 

8x  7  8m  9co  8m 


(3.12) 


with  convection  terms  on  the  left,  conduction  terms  first  on  the  right 

8(f> 

and  all  other  terms  not  involving  8m  lumped  into  one  source  term  "d". 

•  ii  ,ii  d^m  .  ii  dijjp  •  " 

The  terms  and  mr  defined  by  -i —  =  -mT  and  -r —  =  -m 
IE  J  6x  I  dx  E 

represent  the  rates  of  mass  transfer  across  the  I  and  E  boundaries. 

At  a  line  of  symmetry  the  mass-transfer  rate  is  zero  while  at  a  wall 
the  mass  flux  is  derived  from  given  boundary  conditions.  When  the  E 
boundary  corresponds  to  the  "edge"  of  a  boundary  layer,  as  in  the 

#  II  M 

problem  considered  here,  an  expression  for  m.  (=m  )  must  be  derived 

t  b 

to  determine  the  amount  of  mass  entrained  by  the  boundary  layer  at 
its  junction  with  the  free  stream.  The  momentum  equation  (2.2)  in  the 
free  stream  reduces  to 


„  8u  . 

PU  87' 


_  (9P  +  ie  B  ) 
8x  J  z  y; 


(3.13) 


Since  in  the  limit  as  y  approaches  y^, expression  (3 . 13)must 
also  hold  true  just  inside  the  G  boundary,  substitution  back  into  the 
momentum  equation  yields  the  desired  expression  for  the  entrainment 
rate 


1  im 

y-yG 


'  8  /  8u  \  • 
8y  Ueff  8y 

3u/  3y 


(3.14) 
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This  limit  will  exist  and  be  finite  if  the  mixing-length 


concept  is  employed  since  y^^  will  vanish  simultaneously  as  3u/3y 
approaches  zero.  For  example,  substitution  of  Prandtl's  mixing -1  ength 


hypothesis  y  =  p£ 
ment  rate 


3u 


3y 


into  expression  (3.14)  yields  for  the  entrain- 


II 


(3.15) 


It  is  perhaps  worth  noting  at  this  point  that  film-cooling 
mathematical  computational -techniques  empl oying a  simpl e  mixing-length, 
eddy-viscosi ty  hypothesis  have  been  shown  to  be  capable  of  correctly 
describing  the  trends  and  magnitudes  of  important  properties  of  wall- 
jet  and  wall -wake  flows.  Additional  references  include  the  work  of 
Kacker  and  Whitelaw  [65]  who  employed  the  Prandtl -Kol mogorov  model  of 
turbulence  with  a  simple  empirical  expression  for  the  length  scale  to 
obtain  satisfactory  predictions  over  a  wide  range  of  film-cooling  flow 
situations . 


3.2  MHD  Modifications  to  Effective  Viscosity  Hypothesis 

The  effective  viscosity  hypothesis  assumes  that  in  the 
expression  for  the  effective  shear  stress  =  yg^  y^^  is 

composed  of  a  laminar  viscosity  y  plus  an  eddy  viscosity  which, 

•j* 

upon  the  adoption  of  Prandtl's  mixing-length  theory,  may  be  written  as 

e  =  p£^  |  3u/3y  |,  (3.16) 

Various  empirical  relations  have  been  derived  for  the  mixing- 
length  and  many  of  these  can  be  expressed  by  the  general  form 


^Often  Klebanoff's  intermi ttency  factor  T  is  included  in  the  expres¬ 
sion  for  e  for  the  region  away  from  the  wall,  but  it  will  be  assumed 
unity  here. 
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l  =  Dy£f(x,y/y£) 


(3.17) 


where  "D"  is  a  general  damping  function  which  is  zero  at  the  wall  and 


equal  to  unity  at  large  distances  from  the  wall,  "y£"  is  some  charac¬ 
teristic  thickness  of  the  boundary  layer  and  "f"  is  some  function 
which  permits  the  mixing-length  to  vary  with  downstream  and  cross¬ 
stream  di recti ons . 


Several  models  for  have  been  proposed  for  the  region 


near  a  film-cooling  slot  and  attempts  to  account  for  the  effects  of 
relaxation  or  slot-lip  thickness  for  example  have  resulted  in  fairly 
accurate  predictions  by  a  judicious  choice  of  the  function  "f"  and 
appropriate  modifications  of  the  effective  diffusivity  profiles  near 
the  slot  exit.  In  Ref . [66]  abridging  technique  is  employed  as  shown  in 
Fig. 3.1  to  preclude  the  unrealistic  occurence  of  a  vanishing  diffu¬ 
sivity  at  points  of  veloci ty  extremum.  In  Ref. [67]  further  modifications 
are  imposed  on  the  bridged  profile  to  represent  the  effect  of  slot 
lip  thickness  in  the  form  of  an  "additive"  diffusi ti vity  also  shown 
in  F i g . 3 . 1 . 


y 


u 


a )  Velocity  profil  e 


b)  Effective  viscosity  profile 


Fig.  3.1  Bridging  technique  for  effective  viscosity 
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For  Mach  3  and  higher  Mach  number  flows,  Bushnell 
[68]  reports  good  agreement  with  experimental  velocity  profiles  and 
heat  transfer  rates  by  assuming  a  certain  value  for  the  mixing  angle 
a, (Fig. 1. l) and  prescribing  a  mixing-length  proportional  to  the  local 
boundary-layer  thickness  6..  Flow  relaxation  down-stream  of  the 

l  i 

cooling  slot  is  divided  into  three  regions  and  a  maximum  mixing-length 

*max  =  '°8  6i  is  prescribed  where  «i  =  V2’  yG,2  '  yG,l  and  yG  '  yG,l 
respectively.  The  simplified  procedure  is  claimed  to  be  applicable 

to  the  near  slot  (x/y^  <  30)  as  well  as  the  far  slot  (x/y^  >  30) 

regions  even  though  the  boundary-1 ayer  equations  are  not  valid  in  the 

near  slot  region. 

Relaxation  effects  could  be  approximately  accounted  for  in 
Ref.[69]by  assuming  that  the  mixing-length  is  a  function  of  a  velocity 
profile  shape  factor.  According  to  the  work  of  Maise  and  McDonald 
[70]  mixing-length  values  for  equilibrium  incompressible  flows  can 
be  used  for  compressible  equilibrium  flows  with  adiabatic  walls  up  to 
Mach  number  5.  The  function  f(x,y/y^)  was  therefore  chosen  as  a 
function  of  the  incompressible  form  factor  which  is  defined  as  the 
ratio  of  the  displacement  thickness  to  the  momentum  thickness.  Three 
functional  variations  derived  from  experimental  results  were  tried 
with  a  quadratic  variation  giving  the  best  agreement  near  the  end  of 
the  relaxation  region  for  recovery  from  wall  blowing. 

In  the  present  work  the  effective  viscosity,  y  ^  is  assigned 
in  much  the  same  manner  as  was  done  by  Cole  and  Spalding  [3]  with 
further  modifications  introduced  in  an  attempt  to  account  for  the 


. 
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influence  of  the  magnetic  body  force.  The  contribution  to  y  ^  of 
shear-generated  turbulence  is  accomplished  by  adopting  a  mixing-length 
profile  illustrated  in  Fig.  3.2  which  is  appropriate  to  a  velocity 
profile  as  expected  from  a  tangential  injection  slot. 


(a)  Velocity  profile 


(b)  Absolute  tangential-  (c)  Mixing-length 
velocity  gradient 


Fig.  3.2  Effective  viscosity,  mixing-length  model 


3u 

3y  I  mi  n 


determines  characteri s ti c 


In  Fig.  3.2(b)  a  selected  minimum  value 
lengths  (y£)12,  (y£)34  and  (y£)5g.  For  these  far  wall  regions  f(x,y/y£) 
is  assigned  an  empirical  mixing-length  constant  A  and  the  correspondi ng 
mixing-lengths  are  then  £]2  =  DA(y£)12,  £34>  =  DA(y£)34  and  £5g  =  DA(y£)56 
The  value  of  A  is  taken  as  0.075  and  the-  MHD  damping  function  D  will  be 


derived  below.  The  value  of 


Bu 


ay 


mi  n 


is  also  established  from  the  gross 


shape  of  the  velocity  profile  and  equated  to  FR  U^/y^  with  FR  as  another 
empirical  constant  taken  as  0.15  for  best  agreement  with  experiment,  Ref . [3] 
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Equation  (3.16)  used  in  conjunction  with  an  £~y  profile 
such  as  Fig.  3.2(c)  would  result  in  vanishing  e  in  some  regions  of  the 
boundary-1 ayer .  An  unreaslistic  zero  value  for  e  is  precluded  by 
making  an  allowance  for  local  "free  stream"  turbulence  and  postulating 


that 


9u 


9y 


never  falls  below  a  value  typical  of  the  fluctuating  component 


of  the  local  "free  stream"  velocity.  This  minimum  fluctuating  component, 
U1  .  ,  is  assumed  to  be  proportional  to  the  local  velocity,  namely  u'  . 

=  F^u  where  Fu  is  yet  another  empirical  constant  taken  as  0.02. 

Near  the  wall  the  mixing-length  is  prescribed  by  the  usual 
linear  variation  l  =  DKy  with  K  =  0.435.  Here  the  damping  function,  D, 
will  be  a  modified  Van  Driest  function  which  is  derived  in  the  follow¬ 
ing  manner.  To  account  for  the  increasingly  predominant  influence  of 
the  laminar  viscosity,  v,  upon  approaching  a  solid  boundary.  Van  Driest 
[19]  considered  the  velocity  distribution  normal  to  an  oscillating 
semi-infinite  plane  (Stoke's  2nd  problem)  given  by 


u (y )  =  uq  exp 


V(£)y 


COS 


nt -y/(  )y 


(3.18) 


Identifying  u  as  the  fluctuating  velocity  u'  and  taking  the  inverse 
situation  with  the  plane  fixed 

-y/Ai 

u'  =  u  '  (1  -e  ),  (3.19) 


where 


A-|  =  /(2)  v  //(nv) 
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The  quantity  /(nv)  with  units  of  velocity  is  then  replaced  by  the 

"friction  velocity"  /(t  /p  )  and  the  constant  At  is  now  replaced  by 

J  w  w  1  r  J 

A  given  as  A  =  A  v//(x  /p  )  with  A+  chosen  to  agree  with  experi- 


Generalizing  this  damping  effect  to  include  blown  plates  or 
adverse  pressure  gradients  for  which  the  local  shear  stress  can  vary 
markedly  with  y  in  the  region  near  the  wall,  Kays  [71]  among  others, 
argue  that  the  local  shear  stress  is  more  likely  to  effect  the  local 
damping  and  ought  to  be  used  in  the  constant  A  rather  than  the  shear 
stress  at  either  extremity  of  the  boundary  layer.  Hence  the  first 
contribution  to  the  function  D  due  to  viscous  damping  can  be  written  as 


D-j  =  1  -  exp(-y+/(T  +  )/Ao+) 


(3.20) 


where  y+  =  (y/v)/(i/p),  A+  =  (A/v)/(t/p),  and  t+  =  t/t 

U  W 


The  obvious  extension  to  include  MHD  effects  is  to  resolve 
Stoke1 s  2nd  problem  with  the  magnetic  body  force  term  included.  In 
its  simplest  form  the  equation  to  solve  is  ^ 


8u 

3t 


(3.21) 


with  boundary  conditions 

u(t,0)  =  u  cos(nt)  and  u  =  finite  at  y  =  «=,  t  >0 
v  o 


+  Stroke's  first  problem  using  Eq . ( 3 . 21 )  was  first  solved  by  Rossow  [72] 
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The  solution  by  Laplace  transforms  gives 


u  =  uQexp 


"y/(2V  /(m2+n 


cos 


nt-y/(  2vy/(m^+n 


(3.22) 


where 


m  =  aB  /p 
y 


By  reasoning  similar  to  the  one  used  for  viscous  damping 
we  deduce  that  the  second  contribution  to  the  function  D  is  given  by 


°2  =  1 


or 


exp 


-y 


+  /  T 


+ 


V 


A  A 


<mV 


+  1)  + 


«+2  + 
Ao  m  ' 


=  1  -  exp(-y+  */(t"V A+) 


(3.23) 


where  now 


A+  =  A  *//(/( 


+  1)  + 


A +2  + 
A  m 
o 


) 


(3.24) 


and 


m+  =  my/x, 


w 


Finally,  outside  the  viscous  sublayer,  damping  due  to  electro¬ 
magnetic  forces  still  persists  although  viscous  damping  is  negligible. 

-A2m+ 

Resolving  Eq.  (3.21)  without  including  v  results  in  u  =  uQe  from 

which  the  third  contribution  to  the  damping  function  D  is  taken  to  be 

+ 


D3  e 


-A2m 


(3.25) 


Experiments  performed  by  Lykoudis  and  Brouillette  [17,18] 
indicate  that  best  agreement  was  obtained  with  A2  =  700. 

Since  the  v'  fluctuations  are  not  damped  by  the  B  field 
(ie.  they  are  parallel  to  the  magnetic  lines  of  force)  and  since 


=  0 


,  D-j  is  used  to  represent  wall  damping  of  the  v'  fluctua^ 


■ 

tions,  is  used  to  represent  both  wall  and  MHD  damping  of  u'  and 
is  used  to  account  for  MHD  damping  away  from  the  wall.  The  final 
expression  for  D  in  Eq.  (3.17),  therefore,  is  written  as 

D  =  /(D^D^  (3.26) 


One  final  refinement  introduced  in  this  dissertation,  in  an 
attempt  to  account  for  the  shape  alteration  of  the  gross  velocity  pro¬ 
file  by  the  magnetic  body  force,  involves  a  modification  of  the  con¬ 
stant  A  +.  In  OHD  flows,  A  + =  26  is  often  assumed  even  though  it  has 
been  shown  to  be  influenced  appreciably  by  other  factors  such  as  trans¬ 
piration  and/or  pressure  gradients  represented  by  parameters  v  +  and  p+. 
An  empirical  correlation  of  experimental  results  is  given  by  Kays  [71]  a 


A  +  =  A  +  /  ( 5 . 1  5 

o  o,o 


+ 

v  + 


5.86p 
1  +5v  + 


+  1.0) 


(3.27) 


where 


n  +  —  +  ~  r  .  +  +  +  n 
A  A  =26  when  m  =  p  =  v  =  0 
o  o  ,o 


An  analogy  between  wall  suction  and  magnetic  body  force 
influence  on  the  gross  shape  of  the  velocity  profiles  was  first  pointed 
out  by  Lykoudis  [73].  The  simultaneous  presence  of  wall  suction  and 
magnetic  field  was  investigated  further,  by  Luft  and  Rodkiewicz  [40]. 
Equation  (3)  of  this  paper  (which  appears  as  Appendix  A  in  this  work) 
may  be  re-written  in  terms  of  the  variables  defined  here,  namely  m+  and 
y+,  so  that 

-/(m+)  y+ 

u  =  1  -  e  (3.28) 


Comparing  this  with  the  well  known  asymptotic  suction  profile  (Ref.  [74] 
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for  example)  and  defining  v+  =  v  /  (/(t  )/p  )  this  suction  profile 

w  w  c 

may  be  written  as 

i  + 1  + 

-|v  |y 

u  =  1  -  e  (3.29) 


It  is  immediately  apparent  that  /m+  in  (3.28)  plays  an  identical  role 
as  suction  plays  in  the  velocity  profile  given  by  Eq.  (3.29).  Equation 
(10)  of  Ref.  [40]  is  also  seen  to  return  the  asymptotic  suction  profile 
when  the  magnetic  parameter  I  defined  there  is  set  to  zero  and  noting 
that  |v+|y+  e  n  (not  to  be  confused  with  the  definition  of  n  for  film¬ 
cooling  effectiveness  here).  Thus  the  turbulence  parameter  A*  in  the 
functions  D-j  and  is  evaluated  in  accordance  with  Eq . (3.27)  but  with 
v+  replaced  by  /m+  so  that 


5.15 


26 


/m++ 


5 ,86p+ ' 

1  +5/m+ 


+  1 


(3.30) 


CHAPTER  IV 


UNIVERSAL  MHD  MODIFICATIONS  TO  THE  PATANKAR 


AND  SPALDING  NUMERICAL  SOLUTION  PROCEDURE 


4.1  MHD  Wall  Function  Modifications 

The  method  of  Patankar  and  Spalding  involves  a  'Couette 
flow'  analysis  near  a  wall  wherein  the  x-wise  convection  is  neglected, 
by  virtue  of  the  velocity  component  u  being  small.  The  fluxes  of 
mass,  momentum,  and  energy  are  determined  from  the  solution  of  ordi¬ 
nary  differential  equations  which  can  be  solved  once  for  all.  The 
momentum  equation  including  the  MHD  term  may  then  be  written  for  an 
impermeable  wall  as 


or 


y 


where 


o 


and 


In  terms  of  dimensionless  variables  (defined  in  the 


nomenclature)  this  equation  can  be  written  as 


49 


50 


T+  =  1  +  i- 

Tw 


y 


P+Tw 

y„ 


/  ( t  p  )  +  je  B 


wc 


z  y  J 


y 


7  (Vc) 


C  dy+ 


This  suggests  the  definition 


jez  =  — l‘C. 

17  (xw  pc> 


such  that 


y 


+ 


T  =  1  + 


(p+  +  je!)dy+ 


o 


+ 


Since  x  -  y  the  equation  to  be  integrated  is  then 


dy 


+ 


du+  = 
dy+ 


1  + 


y 


o 


(P++je^)dy+ 


(4.1) 


Defining  further  (similar  to  the  Lykoudis  parameter  A  introduced 
in  Ref.  [73]) 


+ 


□  2 

,  a  B  y 
rr  _  a  „+  _  c  y  c 

a  =  —  m  =  - z - 

a  p  x 

C  '  c  w 


and 


be  expressed  as 


B  /  ( x  /p  ) 
y  v  w  c' 


’  Ex  = 

A 


B  /  (x  /p  ) 
y  v  w  cy 


.  4- 

*  Jez 


may 


1+3  L 


Ez+u+-< 


(4.2) 
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For  the  electromagnetic  field  influence  on  the  recovery 

factor  it  turns  out  to  be  more  convenient  to  define  an  electric  field 

parameter  without  reference  to  the  magnetic  field.  Hence  let 

e+  =  E  /  (cru  )/t  and  if  we  further  restrict  the  solutions  to  be 
z  c  c  w 

valid  for  negligible  Hall  effect  in  the  Couette  region  the  current 
density  can  be  expressed  as 

jez  =  m  a  [Ez  +  u  ]  (4.3) 

The  evolution  of  the  parameter  e+,  which  is  related  to  E*  by 
e+  =  /  (m+)Ez  will  be  demonstrated  in  the  following  section  deal¬ 
ing  with  the  integration  of  the  Couette  flow  equations. 

The  global  energy  equation  in  the  Couette  region  may  be 

written  as 


where 


37  =  -  £  {Jh  -  +  jezEz 


,  yeff  9h  _  ,9T  ,  9Te 

n  '  '  Preff  3y  "  K  9y  "  e  9y 


+ 


yp  v  h 
L  s  s  s 
s 


(4.4) 


or,  in  terms  of  the  dimensionless  variables  (defined  in  the  nomen¬ 
clature)  , 


dcp+  _  ^reff  ,  W 
+  +  2 


dy 


y 


(1  +  Pr  ff) 


du 


+2 


Pr 


dy 


+ 


-  2 


eff 


y 


+ 


+ 


y 


jezEzdy 


(4.5) 
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Treating  u+,  Pr  ^  and  c T  as  known  functions  of  y' ,  the  solutions 
to  Eqs .  (4.1)  and  (4.5)  may  be  written  formally  as  functions  of  the 


various  parameters,  namely  u+  =  fu(m+,p+,E*,y+) 

and  y  =  (m+,p+,E^,y+,Pr)  +  M.  >2(m+,p+,E^y+,Pr) 

Specifying  these  solutions  appropriate  to  the  outer  edge  of  the 
Couette  region,  the  rates  of  transfer  of  momentum,  heat,  and  mass 
(in  the  general  case)  across  the  Couette  flow  may  be  expressed  in 
terms  of  a  drag  coefficient  s,  a  Stanton  number  ST,  and  a  recovery 
factor  r.  It  may  be  shown  that 


s  =  f"2  =  t  /(p  u  2) 

u,c  w  VKc  c  ' 


-1 


ST  =  ( f . .  f,  T  J  =  -  J,  / { ( H  . -H  ,)p  u  } 
u,c  cf ,  1 , c x  h  v  ad  w,ad  c  c 


(4.6) 


=  1  ■  <hw,ad-hc)/(uc2/2) 


Finally,  the  parameters  pertinent  to  the  solution  of  the 
equations  above,  may  be  combined  in  such  a  way  that  expressions  for  s, 
St  and  Re  can  be  derived  as  functions  of  the  more  common  dimensionless 
numbers.  These  are:  the  Reynolds  number  of  the  Couette  flow  Re;  a 
pressure  gradient  parameter  F;  the  Hartmann  number  Ha  based  on  the 
thickness  of  the  Couette  layer;  and  a  newly  introduced  electric 
field  parameter,  E.  Their  definitions  are 


■ 
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Re  =  u+Qy+c  =  (puy/y)c 


F  =  P+y+/{u+)  =  (yc/pcuc2)dp/dx 

+  +  +  Ez^ c 

F  =  e  y  /u  e  -  /  (cr  /y  ) 

Jc  c  u  v  c  c' 


(4.7) 


Ha  =  /  (m+)y*  e  Byc  /  (a,/yc) 


Before  the  integration  of  Eqs.  (4.1)  and  (4.5)  is  carried 
out,  it  is  possible  to  eliminate  the  turbulence  parameter  K  by  a 
further  change  of  variables  defined  as  follows: 


y*  =  Ky+ 
u*  =  Ku+ 


<P 


*  = 


Kf 


A*  =  KA 

p*  =  ^ 

p  K 


+ 


m 


★  = 


m_ 

K: 


E*  .  kE; 


je*  =  je*/K 


(4.8a) 
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R*  =  K  Re 
F*  =  F/K2 


Ha  =  Ha 


E*  =  E 


•  *  = 


K 


then  the  following  relations  also  hold 


*  * 

R*  =  u  y 
c  c 


F*  =  (p*y*/u,  ) 


*  'c 


E*  =  e*//m* 


(4.8b) 


Ha  =  / (m*)y* 


E*  =  (e*y*/u*)c 


The  equations  to  be  integrated  then  become 


du*  T 


dy*  + 
7  y 


wi  th 


y1 


T+  =  1  + 


(p*  +  je*)dy* 


(4.10) 


je*  =  m*  [E*+u*] 


(4.11) 


and 


(4.12) 


dj>* 

dy* 


Pr 

eff  W 

+  2K 

y 


(1-Preff)2u* 


du* 

dy* 


Pr 


eff 


y 


+ 


0 


je*E*dy* 
z  z  J  ) 


4.2  Integration  of  the  Couette  Flow  Equations 

Assuming  both  the  pressure  gradient  and  current  density  to 
be  independent  of  y  across  the  Couette  layer,  Eq.  (4.10)  can  be  inte¬ 
grated  to  yield 


+ 


T 


1  +  (p*  +  je*)  y* 


(4.13) 


and  the  curve-fitted  solutions  of  Patankar  and  Spalding  could  be  used 


by  merely  defining  an  effective  pressure  gradient  parameter  such  that 

+ 


eff 


/  +  .  +  \ 
=  (P  +  Jez) 


y. 


/  +\  2 

<uc> 


(4.14) 


This  was  the  approach  adopted  in  Ref.  [7]  although  the  authors  concede 
that  perhaps  je*  should  be  allowed  to  vary  with  y  as  is  evident  from 
expression  (4.11).  In  this  thesis  relations  will  be  presented  which 
do  allow  je*  to  vary  as  prescribed  by  Eq.  (4.11)  although,  as  was 
done  in  the  general  OHD  method  of  Patankar  and  Spalding,  a  constant 
property  Couette  flow  will  be  assumed  such  that  a+  =  1  will  also  be 
imposed . 


Substitution  of  Eq.  (4.11)  into  Eq.  (4.10)  yields 

y* 


(p*  +  m*E*)  y*  + 

m*  u*dy* 

(4.15) 

y*  o 

1  +  P*.^  y*  +  m* 

u*dy* 

(4.16) 

0 

P*4T4r  =  p*  +  m*E* 
eff  z 

(4.17) 

t 


Private  communication. 


which  can  be  interpreted  as  an  effective  pressure  gradient  parameter. 
It  is  evident  from  this  that,  in  the  presence  of  a  magnetic  field 
(m*  f  0),  the  effective  pressure  gradient  is  altered  only  if  an  elec¬ 
tric  field  is  also  present;  that  is,  E*  must  not  be  zero. 

Now,  strictly  speaking,  in  addition  to  Re,  there  are  three 
independent  parameters,  p*,  m*,  and  E*,  which  can  influence  the 
solution  to  Eq.  (4.9).  Fortunately  as  Eq.  (4.16)  shows,  however, 
there  are  effectively  only  two  free  parameters,  namely,  P*^  and  m*, 
since  p*,  m*  and  E*  appear  as  a  group  to  form  a  single  coefficient 
of  y*  in  Eq.  (4.16).  The  equation  for  the  shear  stress,  Eq.  (4.9), 
may  then  be  integrated  and  curve  fit  for  two  special  cases:  one  with 
=  0  and  the  other  with  m*  =  0.  For  m*  =  0  no  electromagnetic 
effects  are  present  and  the  solution  given  by  Patankar  and  Spalding 
for  the  effect  of  pressure  gradient  on  drag,  must  be  recovered.  This 
provided  a  convenient  check  on  the  computer  routine  written  for  the 
drag  coefficient  s.  The  curve-fitted  solution  for  this  case  is 
given  as  [2]. 


4F*R’ 


[12. 82'5  +  R*2-5] 


0.4 


1  .6 


(4.18) 


where  s*  is  the  value  of  s*  as  a  function  of  R*  with  P*ff  =  m*  =  0. 
Then,  if  a  curve-fitted  solution  for  P*^  =  0  could  also  be  found, 
the  resultant  expression  for  s*/s*  as  a  function  of  R*,  p*,  m*  and 
E*  could  be  assumed^  to  be  formed  of  the  product  of  this  solution 
with  Eq.  (4.18)  in  which  F*  would  be  replaced  by  an  effective  pressure 


+This  assumption  must  be  tested  by  comparing  the  results  it  yields 
with  an  exact  integration  which  is  done  in  Table  1. 
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gradient  parameter  F*^  computed  by  replacing  p*  with  P*^,  that  is 

F*ff  =  (p*  +  m*E*)  y*/uj2  (4.19) 

Before  proceeding  with  the  details  of  curve  fitting  the 
solutions  for  P*^  =  0,  the  equation  for  total  enthalpy  with  m*  and 
P*.|r.f  as  parameters  will  now  also  be  written  for  completeness: 


df* 

dy* 


^reff  W 

+  2K 

P 


0-Preff)  2u 


du* 

dy* 


2Pr 


eff 


+ 


y 


e?y 


+ 


u*dy* 


(4.20) 


Furthermore,  substituting  the  results  of  Section  3.2  for  y+ 
prior  to  integration,  the  velocity  equation  (4.9)  becomes 


du*  _  _ 2t^ _ 

dy*  "  2  + 

!+/(!+  4y*  t  D1D2D3) 


(4.21) 


with  Eq.  (4.16)  providing  the  link  for  t+.  Figure  4.1  shows  typical 
results  obtained  from  the  numerical  integration  of  Eq.  (4.21).  The 
curve  of  u*  vs  y*  for  m*  =  0  corresponds  exactly  to  the  curve  labelled 
p*  =  0  in  Ref.  [2]  as  it  should.”*" 

4.3  Curve-Fitting  the  Integrated  Momentum  Equation  -  Shear  Stress 

A  convenient  feature  of  the  method  of  Patankar  and  Spalding  [2] 
for  OHD  flows  is  the  ability  to  condense  the  integrated  Couette  flow 


•J* 

Note  that  in  Ref.  [2],  m*  represents  a  mass  transfer  parameter  and  is 
not  to  be  confused  with  the  magnetic  field  parameter  m*  as  defined 
in  this  thesis. 


0.004 
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Fig.  4.1  Couette  Layer  Dimensionless-Velocity  Profiles  l|^=0 
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equations  into  a  set  of  algebraic  equations.  Such  equations  will 

now  be  derived  in  this  section  for  MHD  flows  with  the  following  results: 

We  require  a  single  equation  for  the  curves  shown  in  Fig.  4.1 

k 

tor  P*^  =  0.  Figure  4.2  is  a  plot  of  s*/s*  against  Ha  for  various 
Reynolds  numbers,  R*.  Analogous  to  the  effect  of  mass  transfer  on 
drag1  illustrated  by  similar  curves  in  Ref.  [2]  it  is  observed  that 
the  horizontal  distance  between  any  two  curves  of  constant  R*  is  approxi¬ 
mately  uniform  for  all  values  of  s*/s*.  This  suggests  that  all  the 
curves  could  be  made  to  collapse  on  to  a  single  curve  if  the  abscissa 
were  multiplied  by  a  function  of  R*. 

In  the  technique  of  curve  fitting  and  parameter  correlation 
(see,  for  example,  Ref.  [75],  [76])  one  usually  strives  to  obtain 
linear  relationships  among  these  parameters.  In  plotting  s*/s* 

k 

against  Ha  on  rectangular  graph  paper,  linear  curves  are,  unfortunately, 
not  obtained.  Referring  back  to  Eq.  (4.15)  it  is  observed  that  the 
magnetic  field  parameter  m*  appears  only  to  the  first  power  and  pre- 

k  *2  2 

vious  formulation  shows  m*  related  to  Ha  through  m*  =  Ha  /y*  .  This 

*2 

suggests  trying  a  plot  of  s*/s*  against  Ha  as  is  done  in  Fig.  4.3. 

It  is  evident  that  the  general  form  of  the  equation  for  these  curves 
for  all  R*  is 


where  a  and  b  must  be  regarded  as  functions  of  R*.  Now,  for  all 

k 

curves  shown,  s*  =  s*  when  Ha  =  0  as  it  should  by  definition;  hence 


^T his  result  further  substantiates  the  analogy  between  wall  suction 
(mass  transfer)  and  MHD  effects  first  proposed  by  Lykoudis,  [73]. 
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Fig. 4.2  Couette  Layer  Drag  Coefficient  vs.  Hartmann  Number,  Pe*^ 
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a  =  1  and  b  =  b(R*) . 

To  evaluate  b  as  a  function  of  R*,  a  straight  line  was 
eventually  obtained  when  plotting  the  slope  b  against  R*  on  log-log 
coordinates  as  is  shown  in  Fig.  4.4.  Evidently,  b  has  the  functional 
relation;  b  =  nR*  m,  and  graphically  we  find  n  =  117.55,  and  m  =  1.345. 
Hence 

-1  q/iq  *? 

(s*/s*)  -I  =  1  -  117.55  R*  JHDHa^ 
v  o 7  emp ,  1 

1  Ha2  (4>22) 

0.0085  R*1 ’345 

If  the  empirical  values  of  s*/sj  computed  from  Eq.  (4.22)  are  now 
plotted  against  the  exact  values  of  s*/s*  taken  from  the  numerical 
integration  of  the  Couette-flow  equations,  a  45  degree  line  is  ex¬ 
pected  if  a  one-to-one  correspondence  is  to  be  obtained.  Figure  4.5 
shows  such  a  plot  and  it  is  observed  that  agreement  is  fairly  good 
over  a  certain  range  of  s*/s*  and  for  Reynolds  numbers  greater  than 
about  100.  For  smaller  values  of  R*  the  curves  deviate  excessively 
from  the  exact  values  and  a  correction  for  small  Reynolds  numbers 
must  be  derived. 

In  order  to  derive  a  suitable  smal 1 -Reynolds  number  correc¬ 
tion,  the  obvious  solution  is  to  plot  the  difference  6  =  (s*/s*) 

o  exact 

*(s*/s*)emp  against  (s*/sJ)emp  for  various  R*  and  hope  for  a 
family  of  straight  lines.  This  procedure  failed  and  a  plot  of  6 

★  9 

versus  Ha  was  attempted  with  equal  lack  of  success.  Finally,  a 
plot  of  6  against  1  -  (s*/sJ)emp  yielded  a  set  of  reasonably  straight 

'Typical  values  of  R*  =  u*y*  are  in  the  range  of  about  10  to  1000 
as  seen  from  Fig.  4.1. 


Fig.  4.4  Curve  Fitting  Parameter '  b'  vs.  Reynolds  Number 
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Fig.  4.5  Couette  Layer  Drag  Coefficient 
Small  Reynolds  No.  Correction 


lines  such  that  the  general  form  of  the  equation  for  6  could  be 
written  as 


6  =  b( R*) 


*0 

Ha^ 


.1.345 


0.0085  R* 

The  new  expression  for  (s*/s*)emp  2  (which  ideally  equals 
(sVs* ’exact5  may  now  be  written  (s*/s!Pemp,l  +  6  or 


(s*/s^emp,2  =  1  + 


*2 

Ha^ 


0.0085  R* 


1.345 


(4.23) 


Evidently  the  function  b(R*)  in  Eq.  (4.23)  must  possess  the  property: 
b(0)  =  1 ,  so  that  the  second  term  remains  finite  (which,  incidentally, 
is  an  obvious  flaw  in  expression  (4.22))  and  also  it  must  somehow  ap¬ 
proach  zero  asymptotically  for  large  values  of  R*.  This  immediately 

_nR*m  -?  ^D?nR*m 

suggests  the  form  b(R*)  =  10  '  =  e  .A  trial  and  error 

solution  of  plotting  b(R*)  on  semi -log  paper  revealed  a  linear  curve 
for  m  =  1.5  from  which  the  constant  n  was  finally  evaluated  as 

_3 

1.31  x  10  .  Hence  the  resulting  empirical  expression  for  the  mag¬ 

netic  field  influence  on  the  drag  coefficient  takes  the  form 


(s*/s*)  n  9  =  1  -  [1 -exp( -3.02  x  10~3R*1,5)] 
v  o'emp,2  l  i  j 

Figure  4.6  illustrates  a  plot  of  (s*/so)exact 


0.0085  R*1 ,345 


(4.24) 


versus  (s*/ s*) „„„  and 
v  o'emp 


agreement  is  deemed  to  be  reasonable. 

It  now  remains  to  deal  with  the  most  general  situation, 
namely  that  for  which  all  the  parameters  R*,  p*,  m*,  and  E*  are  pre¬ 
scribed  independently  and  where,  specifically,  the  pressure  gradient 
is  non-zero  or  the  electric  field  parameter  E*  takes  on  values  other 


(s  /$o  )exact 
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than  those  dictated  by  E*  =  -  p*/m*.+  The  approach  adopted  in  Ref.  [2] 
to  a  similar  situation  implies  that  a  general  formula  be  assumed  by 
forming  the  product  of  Eqs.  (4.24)  and  (4.18)  with  F*  replaced  by 
F*^  defined  in  Eq.  (4.19).  In  Table  1  is  given  a  comparison  of  such 
a  product  solution  with  the  exact  values  for  s*/s*.  The  agreement  is 
seen  to  be  satisfactory  for  a  range  of  representati ve  flow  parameters. 


Table  1 


R* 

Ha 

F* 

eff 

(s*/s*l 

1  7  bcr  exact 

( s*/s*) 

v  o'emp 

%  Deviation 

14.69 

0 

0.0148 

0.890 

0.920 

4.00 

15.94 

0.179 

0.0252 

0.800 

0.850 

6.40 

18.33 

0 

-0.0148 

1  .050 

1  .080 

3.20 

18.76 

0.223 

-0.0178 

1 .020 

1  .070 

5.40 

19.52 

0.050 

-0.0033 

0.970 

1 .020 

5.10 

20.47 

0.224 

0 

0.910 

0.970 

7.60 

51  .44 

0.900 

-0.0138 

0.806 

0.741 

8.17 

86.54 

0.150 

-0.0045 

0.998 

1  .020 

2.46 

101.16 

0.716 

0 

0.883 

0.885 

0.22 

102.97 

0.810 

-0.0275 

1  .080 

1  .010 

6.75 

103.28 

0.510 

-0.0046 

0.960 

0.970 

0.75 

105.11 

1.160 

0 

0.730 

0.710 

2.70 

328.38 

0 

-0.1002 

1  .673 

1  .715 

2.51 

491 .71 

2.680 

0 

0.818 

0.797 

2.62 

503.33 

0.700 

-0.0135 

1.068 

1  .073 

0.43 

903.86 

4.470 

0 

0.763 

0.751 

1.58 

973.34 

0 

0.0528 

0.668 

0.684 

2.42 

1035.05 

6.320 

0 

0.598 

0.586 

2.07 

1047.57 

3.900 

-0.0200 

0.988 

0.956 

3.35 

1107.97 

0 

0.0815 

0.529 

0.532 

0.59 

Expression  (4.24)  is  valid  for  non-zero  pressure  gradients  also  as 
long  as  the  electric  field  E*  satisfies  this  relation. 
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It  is  evident  from  Table  1  that  within  a  deviation  of  plus  or  minus 
10%  the  wall  shear  stress  can  be  calculated  from 


s*/s* 

o  ,emp 


1  - 


4Fe*ffR* 


[12.82-5+R*2-5]0-4  > 


1.6 


x 


1  -  [1 -exp( -3 .02x1 0~3R*^ ‘ 5 ) ] 


*2 

Ha^ 


0.0085R* 


1.345 


(4.25) 


It  should  be  pointed  out  that  in  using  Eq.  (4.25)  the  magni- 

k 

tude  of  the  parameters  such  as  Ha  must  be  limited  so  as  to  avoid  zero 

or  negative  values  of  (s*/sJ)emp  since  this  would  indicate  boundary 

layer  separation.  Calculations  reveal,  for  example,  that  for  F*^  =  0, 
★ 

and  R*  =  25,  Ha  must  be  <_  1  in  order  that  Eq.  (4.25)  be  accurate  to 

k 

within  10%.  If  Ha  =  Byc  /  (cjc/)j  )  is  calculated  for  typical  values 

2  -4 

of  B  =  5  Tesla  (webers/m  ),  yc  =  5  x  10  m,  oQ  =  10  mho/m  and 

_q.  k 

Pc  =  10  kgm/m-sec,  then  Ha  =  0.78.  Hence  Eq.  (4.25)  does  not  pose 
any  physical  limitations  since,  to  the  author's  knowledge,  maximum 
magnetic  fields  which  can  be  generated  by  present  day  super-conducting 
magnets  are  of  the  order  of  7.5  Teslas. 

4.4  Curve  Fitting  the  Integrated  Energy  Equation  -  Recovery  Factor 
In  order  to  formulate  the  effect  on  the  recovery  factor, 
r/r  ,  due  to  the  magnetic  and  electric  fields  in  a  similar  manner  as 
was  done  for  s*/s*,  it  is  evident  that  Eq . (4 . 20 ) ,  with  m*  and  P*^  as 
parameters,  must  be  re-written  since,  as  expression  (4.27)  shows,  E* 
cannot  be  assigned  different  values  via  P*^  with  m*  =  0.  Further¬ 
more,  with  m*  appearing  in  the  denominator,  integration  of  Eq.  (4.20) 
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cannot  be  performed  with  m*  =  0. 

Recalling  Eq.  (4.12),  the  integral  jezE*dy  can  be  re¬ 
written  in  the  following  way.  First,  the  integrand  may  be  written 
as 


je*E*  = 


o  y  E 

c  c  z 


+ 


°cwcEzBu 


T 


W 


T 


W 


Next,  by  introducing  a  previously  defined  electric  field  parameter 


E* 


E  y  a 

(— )  = 
u  y 


e*y* 

u* 


we  obtain 


je*E* 


(m*)e*u* 


(4.26) 


The  simultaneous  equations  to  be  integrated  then  become  Eq.  (4.21) 
with  t  given  by 


x 


+ 


1  +  ( p*  +  e*  /  m*)y*  +  m* 


u*dy* 


(4.27) 


and  the  energy  equation  by 


d(f>*  _  ^reff  W 

dy*  +  2K 

y 


(1  -  Preff)  2u*  ££ 


dy' 


2Pr 


[e*2y*  +  e*  /  rri* 


+ 


u*dy*] 


(4.28) 


Figure  4.7  shows  integration  results  with  p*  =  m*  =  0  when  r/rQ  is 
plotted  against  E*2.  Curve  fitting  these  results  gives  the  following 
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empirical  expressions: 

r/rQ  =  1  +  exp(-0.01875R*)(E*)2  R*  <55 

r/rQ  =  1  +  24(R*)_1 •04(e*)2  R*  >  55 


(4.29) 


In  Table  2  are  given  the  results  obtained  from  expressions  (4.29) 
and  compared  with  the  exact  values. 


Table  2 


R* 

+E* 

^  f//ro^exact 

(r/r  ) 
v  o'emp 

% 

Deviation 

11 

3.33 

10.00 

10.03 

0.3 

53 

5.69 

12.57 

12.99 

3.3 

102 

2.83 

2.58 

2.56 

0.8 

200 

13.49 

18.91 

18.67 

1.3 

527 

0.93 

1  .04 

1 .03 

1  .0 

1060 

47.75 

41  .16 

40.07 

2.7 

Integration  results  with  non-zero  pressure 

gradient  (p*  /  0) 

and  magnetic  field  (m*  /  0)  revealed  r/rQ  to  be  a  weak  function  of  p* 

+  * 

and  not  a  direct  function  of  Ha  when  E*  =  0.  For  small  non-zero 

values  of 

E* 

■k 

the  i nfl uence  of  Ha  on 

r/rQ  appears  to 

be 

negligible  so 

that  Eq. 

(4.19)  can 

be  considered  to 

account  for  the 

primary  electro- 

magneti c 

effects  on 

the  recovery  factor..  For  large 

electric  fields. 

however,  such  as  those  encountered  in  MHD  accel erators ,  a  further 
modification  to  these  expressions  would  be  necessary. 


4- 

‘Note  however  that  even  in  this  case,  r/r  will  be  an  indirect  func¬ 
tion  of  Ha  since  this  parameter  influences  the  wall  shear  stress 
which  in  turn  affects  the  viscous  dissipation. 
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E* 

Fig.  4.7Couette  Layer  Recovery  Factor  vs.  Electrical  Field  Parameter  ( m  =  0 ) 


CHAPTER  V 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 
5.1  General  Discussion 

In  order  to  assess  the  fundamental  influences  of  MHD  body 
forces  and  Joulean  dissipation  on  the  adiabatic  film-cooling  effec¬ 
tiveness,  several  computations  were  performed  by  assuming  thermal 
and  ionizational  equilibrium  as  discussed  in  Section  5.3.  This  assumes 
that  the  electron  gas  is  at  the  same  temperature  as  the  plasma  and 
that  the  mass  electron  number  density  can  be  calculated  from  the 
Saha  relation,  Eq.  (2.7).  These  results  will  also  be  required  for 
assessing  the  effects  of  thermal  non-equilibrium  and  finite  rates  as 
discussed  in  Section  5.4. 

The  MHD  effects  on  the  adiabatic  film-cooling  effectiveness, 
q,  are  clearly  shown  in  the  figures  discussed  in  the  following 
Sections.  However,  to  research  and  understand  the  causes  for  the 
observed  effects,  it  becomes  necessary  to  examine  a  more  intimate 
picture  of  the  phenomena;  examination  of  the  velocity  and  temperature 
profiles  across  the  boundary  layer  helps  to  provide  this  insight. 

These  profiles  as  well  as  corresponding  electric  current  and  plasma 
conductivity  profiles  are  examined  in  detail  in  the  following 
Sections . 

Previous  research  into  OHD  film-cooling  has  shown  that  the 
process  depends  on  many  factors.  Examples  of  these  would  be  the 
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geometry  of  the  film-cooling  injection  slot,  the  properties  of  the 
primary  and  secondary  (coolant)  flows,  the  rate  of  coolant  blowing 
represented  by  a  blowing  parameter  C^p^U^/pgU^,  ratios  such  as  those 
of  temperatures ,  T^/T^  (or  densities)  and  velocities,  U^/U^,  magnitude 
of  free  stream  turbulence  intensity,  magnitudes  of  favourable  or  adverse 
pressure  gradients,  chemical  composition  of  coolant  fluid  etc.  Clearly, 
a  reliable  theoretical  prediction  procedure  can  be  an  expedient 
computation  tool  to  assess  the  effects  of  many  of  these  parameters, 
without  resorting  to  costly  and  time  consuming  experimental  measurements, 
by  simply  changing  a  few  program  statements  or  data  cards.  Any  theoretical 
computational  method  must,  however,  agree  to  within  acceptable  limits 
with  experimental  measurements  whenever  these  are  available. 

Section  5.2  assesses  the  reliability  of  the  present  numerical 
procedure  by  comparing  its  predictions  with  experimental  results  as 
well  as  with  predictions  obtained  with  a  similar  procedure  when 
applied  to  OHD  film-cooling  problems. 

For  the  numerical  solution  techniques,  it  also  becomes 
important  to  investigate  the  sensitivity  of  predictions  to  the  initial 
conditions.  These  are  assumed  to  exist  at  the  beginning  of  the 
numerical  marching  procedure  and,  in  this  case,  represent  assumed 
conditions  at  the  film-cooling  slot  exit.  If  results  prove  to  be 
strongly  dependent  on  the  prescribed  initial  conditions,  then  it 
becomes  necessary  to  solve  the  fully  developed  MHD  channel  flow 
upstream  of  the  slot  and  then  adopt  these  solutions  to  obtain  the 
initial  conditions  at  the  slot  exit.  Section  5.3.4.  includes  an 
analysis  of  the  influence  of  initial  magnetic  field  variation  on  the 
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boundary  layer  velocity  and  temperature  profiles  and  hence  on  the 
film-cooling  effectiveness. 

Section  5.4  considers  the  influence  of  thermal  non-equili¬ 
brium  and  finite-rate  ionization  on  the  MHD  film-cooling  effectiveness. 
Core  flow  or  free-stream  conditions  are  chosen  such  that  little  or  no 
electron  thermal  non-equilibrium  would  exist.  This  insures  that  any 
differences  from  the  equilibrium  solutions  will  be  boundary  layer 
effects  and  not  due  to  free  stream  effects. 

Furthermore,  in  order  to  preclude  complications  for  the  non¬ 
equilibrium  analysis  resulting  from  Hall  currents,  je  ,  the  assumption 

X 

of  infinitely  segmented  electrodes  will  be  made  such  that  je  =0 

X 

everywhere.  To  assess  Hall  current  effects  in  the  boundary  layer, 
this  restriction  is  later  replaced  with  the  assumption  of  free-stream 
Hal  1 -neutral ized  flow,  je  r=0  and  finally  the  condition  E  =0  is 
imposed  which  represents  an  MHD  channel  for  which  the  Hall  potential 
is  short-ci rcui ted. 

5.2  Evaluation  of  Computational  Procedure  for  OHD  Film-Cooling  Flows 
For  the  purposes  of  testing  the  lengthy  computer  program 
written  for  this  thesis,  the  film-cooling  effectiveness  predictions 
derived  from  it  for  ordinary  hydrodynamic  flows  (OHD)  were  compared 
both  with  similar  calculation  techniques  [3],  and  with  experi¬ 
mental  data,  [64]. 

The  present  computer  method  was  compared  with  that  of 
Cole  and  Spalding  [3]  for  which  their  "standard  case"  results  are 

shown  for  comparison  in  Fig.  5.1.  The  initial  flow  parameters  are 
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Fig.  5.1  Validation  of  Computer  Program  for  OHD  Film-Cooling  Effectiveness 
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detailed  in  a  tabular  form  in  Fig.  5.1.  It  can  be  seen  that  both 
computer  predictions  are  virtually  identical  from  which  it  may  be 
concluded  that  the  program  is  functioning  properly. 

Additional  validation  and  credence  to  the  numerical  tech¬ 
nique,  in  general,  is  provided  by  favourable  comparison  with  the 
experimental  measurements  of  Burns  and  Stollery  [64]  for  different 
combinations  of  main  and  coolant  gases.  Figure  5.1  displays  results 
for  helium  and  arcton-12  coolant  injection  at  different  coolant- 
main-gas  velocity  and  density  ratios.  It  is  concluded  that  agreement 
between  prediction  and  experiment  is  sati sfactory . 

5.3  Thermal  and  Ionizational  Equilibrium  MHD  Film-Cooling 

Many  non  film-cooled  MHD  flows  have  been  studied  under  the 
assumption  of  thermal  equilibrium  between  the  electron  gas  and 
carrier  gas  temperature  with  the  electron  number  density,  ng  (or  mass 
fraction  C  ),  calculated  from  the  Saha  equilibrium  relation.  This 
assumption  is  generally  valid  for  open-cycle  combustion  MHD  generators. 
Equilibrium  conditions  have  also  been  assumed  in  this  Section  so  that 
the  influence  of  electromagnetic  body  forces  and  Joule  dissipation  on 
the  film-cooling  effectiveness  can  be  assessed  in  the  light  of  general 
MHD  effects  on  plasma  channel  flows. 

5.3.1  MHD  Influence  on  Film-cooled  Boundary  Layer  Velocity  Profiles 
and  Film-Cooling  Effectiveness 

The  close  agreement  achieved  both  with  the  prediction  method 
of  Cole  and  Spalding  [3]  and  with  experiment  as  shown  in  Fig.  5.1  is, 
in  part,  due  to  having  chosen  identical  initial  velocity  profiles  at  the 
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film-cooling  slot  exit  in  the  prediction  method.  This  initial  profile 
was  also  chosen  for  the  developing  MHD  flow  field  and  is  shown  plotted 
in  Fig.  5.2a  as  represented  by  x/y^=0. 

Many  computer  runs  were  performed  with  different  combinations 
of  the  parameters  available,  such  as  B,  Ur/Ur  ,  Q,  etc.  A  typical  case 
is  represented  in  Fig.  5.2  where  two  identical  runs  were  performed 
except  that  one  is  for  OHD  flow  (B=0)  and  the  other  is  for  MHD  flow 
(ie.  B=2).  The  run  conditions  are  detailed  on  Fig.  5.2a,  which  compares 
developing  velocity  profiles  in  MHD  and  OHD  flows  in  the  near-slot  region. 
Similar  velocity  profiles  for  larger  down-stream  distances  are  plotted 
in  Fig.  5.2b. 

The  developing  velocity  profiles  for  OHD  film-cooling  are 
typical  of  what  is  obtained  by  other  procedures  and  experimental 
measurements.  Spalding  [51]  makes  note  of  developing  velocity  profiles 
becoming  gradually  fuller  "...  but,  that  even  at  100  slot-heights 
downstream,  the  profiles  have  not  lost  all  trace  of  their  original 
depression."  A  similar  trend  is  observed  with  the  OHD  profiles 
shown  in  Fig.  5.2  where  the  velocity  profiles  appear  not  to  be  fully 
developed  even  at  larger  downstream  distances  for  these  run  conditions. 

It  can  be  seen  that  the  OHD  velocity  profiles  lose  their 
maxima  relatively  quickly  compared  to  those  for  MHD  flow  which  are 
seen  to  develop  the  characteristics  of  a  sustained  wall  jet  flow. 

The  MHD  velocity  profiles  are  showing  a  relative  acceleration  of  the 
boundary  layer  fluid  with  that  of  the  free  stream  (core  flow),  since 
the  j”  X  ?  retarding  force  is  less  intense  in  the  boundary  layer  where 
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Fig. 5.2a  MHD  Influence  on  Film- Cooled  Boundary  Layer  Velocity  Profiles-Near  Slot  Region 
Equilibrium  Case 
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Fig. 5. 2b  MHD  Influence  on  Film-Cooled  Boundary  Layer  Velocity  Profiles  -  Far  Slot  Region 
Equilibrium  Case. 
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lower  temperatures  result  in  lower  plasma  electrical  conductivity. 

The  correspondi ng  electrical  current  and  conductivity  profiles  at  the 
same  downstream  distances  are  shown  in  Fig.  5.5  and  Fig.  5.7,  respec- 
ti vely . 

The  MHD  velocity  profiles  imply  that  as  the  magnetic  field 
strength,  B,  is  increased,  more  mass  flow  in  the  channel  gets 
confined  in  the  boundary  layers  of  the  insulator  walls.  The  sus¬ 
taining  of  a  wall  jet  type  flow  by  the  magnetic  body  forces  would  be 
expected  to  have  a  beneficial  influence  on  the  film-cooling  effec¬ 
tiveness,  at  least  for  smaller  values  of  B. 

With  reference  to  Fig.  5.2,  the  wall  jet  can  be  imagined  as 
originating  near  the  mid  point  of  the  upper  slot  lip  and  emerging 
into  a  quiescent  fluid.  The  jet  flow  would  act  as  a  barrier  to 
mixing  of  the  primary  (core)  and  secondary  (coolant)  gases.  At 

larger  values  of  B,  however,  the  higher  j  X  B  forces  in  the  free 

stream  would  direct  most  of  the  channel  mass  flow  into  the  boundary 
layers  of  the  B  walls.  The  hotter  core  flow  gas  would  therefore  be 
forced  into  the  boundary  layer  in  sufficient  quantity  to  rapidly 
heat  up  the  coolant,  resulting  in  a  decrease  in  effectiveness.  It 
would  also  be  expected,  by  the  observed  large  difference  in  velocity 
profiles  between  OHD  and  MHD  flow,  that  the  difference  in  n  would  also 

be  relatively  large  although  the  difference  in  n  from  one  non-zero 

value  of  B  to  another  would  be  expected  less  pronounced. 

The  film-cooling  effectiveness  is  shown  on  Fig.  5.3  and 
plotted  for  different  values  of  magnetic  field,  B.  It  can  be  seen 
that  the  expected  behaviour  of  n  from  studying  the  velocity  profiles  of 
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Fig. 5. 3  MHD  Influence  on  Adiabatic  Film  -  Cooling  Effectiveness  —  Equilibrium  Case. 
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Fig.  5.2,  is  borne  out  by  the  observed  behaviour  of  n  which  indeed 
supports  the  predicted  observations.  All  runs  were  performed  with 
the  flow  parameters  as  shown,  and  only  the  magnetic  field,  B,  was 
varied.  Curves  plotted  are  for  B=0,  1,  2,  and  7. 

It  is  seen  that  the  increase  in  n  from  B=0  to  B=1  is  rela¬ 
tively  much  greater  than  that  for  B=1  to  B=2.  In  fact,  the  predicted 
decrease  in  n  for  larger  values  of  B  is  borne  out  by  the  curves  for 
B=2  and  B=7.  The  latter  actually  shows  that  the  magnetic  field  re¬ 
tarding  force  for  B=7  is  sufficiently  strong  enough  to  cause  flow 
separation  at  approximately  55  slot  heights  downstream  from  the  slot 
exit. 

The  optimum  magnetic  field  intensity  resulting  in  the 
largest  increase  in  film-cooling  effectiveness  is  seen  from  Fig.  5.3 
to  occur  for  values  of  B  between  1  and  2.  By  way  of  comparison,  the 
OHD  film-cooling  effectiveness  is  reduced  to  50 t  at  120  slot  heights 
downstream  whereas  a  similar  reduction  by  50%  does  not  occur  until 
approximately  250  slot  heights  downstream  have  been  reached  in  the 
MHD  case. 

The  influence  of  the  magnetic  field  damping  of  the  turbulent 
velocity  fluctuations  on  the  film-cooling  effectiveness  is  also  ill¬ 
ustrated  in  Fig.  5.3.  Because  any  damping  of  turbulent  fluctuations 
would  reduce  the  rate  of  mixing  of  the  primary  and  secondary  streams, 
the  film-cooling  effectiveness  would  be  expected  to  increase  and  this 
is  verified  in  Fig.  5.3.  The  average  increase  in  n  for  a  magnetic 
field  parameter  B=2  is  seen  to  be  15%  for  the  run  conditions  shown. 
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The  plot  of  the  velocity  profiles  in  Fig.  5.2b  also  shows  the 

locus  of  velocity  maxima  in  the  MHD  boundary  layer  flow.  The  gradual 

decreasing  trend  of  this  curve  can  be  attributed  to  the  increase  of 

electron  current  density,  jez,  at  closer  distances  from  the  walls. 

These  currents,  in  turn,  result  from  higher  plasma  conductivities 

which  are  strong  functions  of  temperature  and  electron  number  density 

which,  as  seen  from  Fig.  5.4,  are  increasing  at  points  closer  to  the 

wall  as  the  downstream  flow  develops.  (The  natural  variation  of  Tg 

with  temperature  (or  degree  of  ionization)  behaves  according  to  the 

following  pattern:  For  degrees  of  ionization  less  than  0.1%,  Tg  is 

an  exponential  function  of  T,  rising  rapidly;  as  the  ionization 

exceeds  this  value,  Tg  ^  T  ).  The  equilibrium  electron  mass-fraction 

profiles,  C*  /  C  ,  are  plotted  in  Fig.  5.6  for  similar  run  conditions 
e  eG 

and  with  B=0  and  B=2.  All  computational  runs  presented  in  this  work 
include  electromagnetic  field  damping  of  turbulent  fluctuations. 

5.3.2  MHD  Influence  on  Film-Cooled  Boundary  Layer  Temperature 
Profiles  and  Film-Cooling  Effectiveness 

The  characteri sties  of  a  wall  jet  flow  for  the  MHD  velocity 
profiles  shown  in  Fig.  5.2  would  result  in  the  upper  half  of  the  jet 
drawing  into  it  the  hotter  fluid  of  the  core  flow  while  the  lower 
half  of  the  jet  would  do  the  same  for  the  coolant  fluid. 

As  discussed  in  the  previous  Section  5.3.1,  the  effect  of 
this  jet  in  keeping  two  fluid  streams  "separated"  for  a  considerable 
distance  downstream,  results  in  a  reduced  heat  transfer  between  them 
since  the  rate  of  mixing  is  reduced  and  velocity  fluctuations  are 
damped.  This  would  result  in  lower  boundary  layer  temperatures  near 
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Fig.  5.4  MHD  Influence  on  Film-Cooled  Boundary  Layer  Temperature  Profiles- Equilibrium  Case 
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the  wall  and  hence  an  increased  effectiveness.  Both  these  observa¬ 
tions  are  substantiated  from  Figs.  5.3,  5.4  where  the  expected 
trends  are  easily  seen.  With  the  run  conditions  used  for  Fig.  5.4, 
the  reduction  in  boundary-1 ayer  temperatures  near  the  wall,  when 
going  from  OHD  flow  to  MHD  film-cooling  (B=2  in  this  case),  is 
approximately  10%  at  downstream  distances  of  interest.  The  correspond¬ 
ing  increase  in  film-cooling  effectiveness  is  seen  from  Fig.  5.3,  and, 
for  a  particular  downstream  position  of  x/yc=200,  amounts  to  approx¬ 
imately  55%.  The  temperature  profile  given  for  twenty  slot-heights 
in  Fig.  5.4  shows  the  effect  of  convective  heat  transfer  resulting 
from  the  core  flow  mass  entrained  into  the  mixing  region  as  depicted 
by  the  velocity  profile  depression  for  x/y^=20  and,  at  the  same 
normal  distance  from  the  wall,  as  seen  by  comparing  Figs.  5.2  and  5.4. 
The  core  flow  temperature  difference  for  a  zero  and  non-zero  (B=2) 
magnetic  field  are  not  significant,  so  that  the  normalized  temperature 
curves  in  Fig.  5.4  can  be  compared  directly  as  actual  temperature 
differences.  For  example,  at  x/y^,=150,  the  free  stream  temperature, 

Tq,  for  the  OHD  flow,  is  2995°K  and  for  B=2  this  value  is  2946°K  -  a 
nominal  2%  lower  value. 

Because  of  some  primary  and  secondary  flow  mixing  and  resulting 
convection  heat  transfer  as  well  as  heat  transfer  through  conduction, 
the  temperature  at  increasing  distances  from  the  slot  increases  at 
closer  points  to  the  wall.  This  trend  is  clearly  evident  on  Fig.  5.4 
which  results  in  part  from  the  decreasing  locus  of  the  velocity  maxima 
as  discussed  in  Fig.  5.2b. 
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Fig. 5.5  Electron  -  Current  Density  Boundary  - 
Layer  Profiles,  Equilibrium  Case 
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Fig.  5.6  Equilibrium  Electron  Mass  Fraction  Boundary-Layer  Profiles 
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Fig. 5. 7  Equilibrium  Plasma  Electrical  Conductivity  Boundary- Layer  Profi les 


The  temperature  boundary  layer  profiles  observed  for  the 
run  conditions  of  Fig  5.4  do  not  display  any  sign  of  Joule  dissipa¬ 
tive  heating,  even  at  350  slot  heights  downstream.  The  profiles 
show  that  the  temperature  in  the  film-cooled  boundary  layer  is 
everywhere  less  than  that  of  the  free  stream.  Joule  dissipation 
effects  are,  of  course,  more  pronounced  in  electrode  boundary  layers 
where  high  electric  currents  prevail.  Viscous  dissipative  heating 
is  seen  to  be  insignificant  as  well  for  these  run  conditions. 

Actual  MHD  generator  channel  3-walls  will,  of  course,  not 
be  completely  adiabatic  especially  when  high  temperature  gas  flows  are 
encountered.  The  heat  transfer  experienced  in  actual  channel  walls 
would  reduce  boundary  layer  temperatures  so  that  Fig  5.4  represents 
the  maximum  possible  temperatures  attainable  in  the  boundary  layer. 

As  mentioned  earlier,  this  is  the  most  important  information  required 
by  the  design  engineer  for  the  generator  design  although  the  heat 
transfer  rate  would  have  to  be  known  to  predict  MHD  generator  per¬ 
formance  and  efficiency. 

In  comparing  Fig  5.2,  5.3,  and  5.4  together  it  is  apparent 
that,  for  equilibrium  plasma  flow,  the  major  cause  for  the  observed 
increase  in  the  film-cooling  effectiveness  over  that  in  OHD  flows, 
is  the  significant  alteration  of  the  boundary-1 ayer  velocity  profiles 
with  appreciable  contributions  to  higher  n  resulting  from  MHD  turbulen 
damping  and  reduced  boundary  layer  temperatures  at  the  wall. 
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5.3.3  Influence  of  Variation  in  Flow-Parameters  on  MHD  Film-Cooling 
Effectiveness  at  Constant  Magnetic  Field 

As  mentioned  in  Section  1.2.3,  extensive  research  both 
experimental  and  theoretical,  has  been  undertaken  for  OHD  film¬ 
cooling  to  investigate  the  influences  of  different  flow  parameters 
on  the  film-cooling  effectiveness.  Some  of  the  studies  relevant 
to  the  tangential -injection  slot  geometry  (chosen  for  the  present 
investigation)  are  briefly  discussed  here.  Whitelaw  [77]  presented 
measurements  of  the  impervious  wall  effectiveness  of  a  two  dimen¬ 
sional  wall  jet.  The  measurements  were  obtained  using  slot  heights 
of  0.5,  0.25,  0.132  and  0.074  inches,  several  si ot-to-free-stream 
velocity  ratios  Uq/Uq  ranging  from  0.288  to  2.66,  and  coolant 
Reynolds  numbers,  745  £  Re^  £  44,000.  The  measurements  confirmed 
results  from  other  investigators  that  the  optimum  effectiveness  is 
obtained  when  Uq/Uqq  is  close  to  unity  but  that  the  magnitude  of  r) 
is  dependent  upon  the  slot  height  to  a  considerable  extent. 

Eckert  and  Birkebak  [78]  studied  the  effect  of  slot  geometry 
on  film-cooling  in  a  system  which  used  air  in  the  primary  and 
secondary  flows.  Experiments  were  made  with  a  tangential -injection 
slot  and  velocity  and  temperature  boundary-layer  profiles  were  mea¬ 
sured  from  which  the  film-cooling  effectiveness  was  calculated. 

Wider  slots  gave  better  performance  in  the  sense  that  less  mixing 
occurred  at  the  slot  exit  resulting  in  somewhat  higher  values  of  n- 
Results  when  compared  to  measurements  of  Seban  and  Back  [79,80] 
indicated  that  the  optimum  slot  geometry  depends  on  the  ratio  of  the 
mass  velocities  in  the  main  (primary)  stream  and  in  the  coolant 
(secondary)  stream  and  also  on  its  range  of  variation. 


Cole  and  Spalding  [3]  applied  their  finite-difference 
numerical  procedure  to  investigate  the  influence  of  the  initial 
velocity  profile  shape  and  slot  lip  thickness,  t,  ,  on  r\  and  made 
comparisons  with  the  experimental  measurements  of  others.  The  shape 
of  the  initial  velocity  profile  was  found  to  have  a  weak  influence  on 
the  adiabatic  wal  1  effectiveness  although  a  more  realistic  initial 
velocity  profile,  also  prescribed  in  the  present  work,  gave  better 
agreement  with  experiment  for  the  near  slot  region  over  the  range 
of  velocity  ratios  computed,  0.5  <  LL/IL  <  1.7.  The  effect  of  t, 

—  L  uO  —  L 

was  not  completely  resolved  and  further  refinement  was  planned  for 
improved  comparison  with  experimental  measurements. 

Measurements  made  by  Kacker  and  Whitelaw  [81]  related  to 
the  influence  of  initial  velocity  profiles  at  the  tangetial -injection 
slot,  supported  conclusions  that  such  film-cooled  boundary  layers 
have  short  lived  memories.  The  effect  on  the  impervious  wall  effect¬ 
iveness  of  the  thickness  of  the  initial  boundary  layer  on  the  upper 
lip  of  the  slot  was  shown  to  be  not  greater  than  5  per  cent.  The 
range  of  velocity  ratios  considered  was  0.3  £  U^/U^  — 

The  same  authors  [65]  also  used  a  numerical  procedure  for  solving 
the  steady,  two-dimensional,  constant  property  form  of  the  Navier- 
Stokes  equations  to  obtain  predictions  of  mean  and  fluctuating  pro¬ 
perties  downstream  of  a  tangential -injection  film-cooling  slot. 

Among  the  parameters  studied,  which  influence  film-cooling  effective¬ 
ness,  were  the  velocity  ratio,  Reynolds  number,  upper  lip  boundary- 
layer  thickness,  shape  of  the  slot  velocity  profiles,  slot  and  free- 
stream  turbulence  intensities,  and  the  slot  lip  thickness. 
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In  order  to  assess  the  correspond! ng  influence  of  some  of 
these  parameters  on  the  MHD  film-cooling  effectiveness,  several  runs 
were  performed  where  the  magnetic  field  was  kept  constant  (B=2  in 
this  case)  and  other  parameters  were  varied.  Figure  5.8  shows  the 
film-cooling  effectiveness  at  different  magnitudes  of  primary  and 
secondary  flow  velocities  and  temperatures .  Two  computations  are 
also  shown  at  similar  temperatures  and  velocities  but  with  B=0  to 
illustrate  the  relative  influence  of  B  compared  with  the  effects  of 
other  flow  parameters. 

It  can  be  seen  that  the  influence  of  velocity  and  tempera¬ 
ture  on  the  OHD  film-cooling  effectiveness  (curves  6  and  7)  is  much 
more  pronounced  than  the  correspond!' ng  influence  of  the  same  parameters 
on  the  MHD  film-cooling  effectiveness  (curves  1  and  4).  In  other 
words,  the  existence  of  electromagnetic  body  forces  suppresses  the 
magnitude  of  the  influence  that  these  same  parameters  exert  on  n  in 
OHD  film-cooling  flows. 

It  is  observed  furthermore,  that  the  magnitude  of  change  in 
n  from  OHD  to  MHD  film-cooling  effectiveness  depends  significantly  on 
the  magnitudes  of  the  initial  velocities  and  temperatures  of  the 
primary  and  secondary  flows  (curves  1  and  7  and  4  and  6). 

The  difference  in  magnitude  as  well  as  in  direction  depends 
also  on  the  downstream  position  of  the  region  in  question.  For 
example,  at  200  slot  heights  the  increase  in  n  from  OHD  to  MHD  flow 
was  only  about  18%  for  Uqo=250,  U^=125,  TgQ=2500,  T^=l 250  (curves  1 
and  7)  while  the  correspond! ng  increase  in  n  was  approximately  250% 
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Fig. 5.8  Influence  onMHD  Film- Cooling  Effectiveness  for  Constant  Band  Variable 
Core  Flow  Parameters-  Equilibrium  Case 


'f°r  Uqo=580,  Uq=290,  Tgo=3000,  T^=1500  (curves  4  and  6)!  For  both 
run  comparisons  it  is  observed  that  in  the  near-slot  region  (small 
x/Yq)  the  film-cooling  effectiveness  is  reduced  for  MHD  flows.  The 
cross  over  point  occurs  at  x/y^  -  80  for  the  lower  initial  velocities 
and  temperatures  (curves  1  and  7)  while  a  similar  change  occurs  at 
x/y^  -  60  for  the  higher  initial  values  (curves  4  and  6). 

The  results  shown  on  Fig.  5.8  also  indicate  that  for  a  specifi 
value  of  magnetic  field  (B=2  in  this  case),  the  film-cooling  effect¬ 
iveness  remains  a  fairly  constant  function  of  x/y^  after  certain 
values  of  initial  velocity  and  temperature  have  been  reached.  For 
example,  at  x/y^=200  the  increase  in  MHD  film-cooling  effectiveness  was 
approximately  19%  in  changing  from  Uqq=250,  U^= 125,  TgQ=2500,  T^=l 250 
to  Ugo=350,  11^=175,  Tgo=3000,  T^=1500  (curves  1  and  3).  Any  further 
increase  in  these  parameters  resulted  in  a  negligible  change  of 
film-cooling  effectiveness  (curves  2,  3,  and  4). 

A  final  observation  made  in  Fig.  5.8  concerns  curve  3  with 

Up  =1000,  Ur=500,  Tp  =3500,  Tp=2000,  and  B=2.  The  effectiveness  is 
Go  C  Go  C 

seen  to  drop  off  rapidly  to  zero  at  approximately  200  slot  heights 
downstream.  This  would  indicate  that  the  flow  has  become  choked  at 
x/y^  -  200  which  is  a  consequence  of  the  physical  channel  size  and 
slot  geometry  chosen  in  the  present  work  and  the  fact  that  the 
initial  Mach  number  was  already  close  to  unity  (MGq=0.86). 


5.3.4  Influence  of  Initial  Magnetic  Field,  Bq(x),  Variation  on 
Boundary  Layer  Profiles  and  Film-Cooling  Effectiveness 

Because  of  computational  problems  arising  in  the  finite- 

difference  technique,  Sherman,  et  al .  [8]  were  forced  to  allow  a 
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gradual  but  rapid  increase  in  axial  variation  of  the  assumed  initial 

magnetic-field  B  (x)  .  Similar  difficulties  arose  also  in  the  present 

investigation  where  computational  problems  were  encountered  in  using 

a  step  function  for  B  (x)  at  x=0. 

o 

It  has  been  shown  in  OHD  film-cooling  research,  for  example 
Refs.  [3]  and  [65],  that  the  shape  of  the  initial  boundary  layer  pro¬ 
files  have  relatively  small  effects  on  p  at  downstream  distances  of 
main  interest.  Computations  performed  in  the  present  work  indicate  that 
similar  conclusions  can  be  stated  for  the  MHD  film-cooling  effective¬ 
ness.  The  functional  relationships  for  B  (x),  assumed  as  initial 
condition  for  the  MHD  finite-difference  procedure,  however,  have 
shown  to  be  of  more  significance.  Several  runs  were  therefore 
performed  with  a  linearly  increasing  function,  B  (x),  whose  slope  was 
varied  in  a  gradual  approach  to  a  step  function.  Figure  5.9  compares 
the  developing  velocity  profiles  for  the  identical  run  conditions 

shown  except  that  (x/yr)  =1  and  10,  respecti vely ,  in  the  relation 

c  ma  x 

Bo=Bmax  at  x/yC=(x/yC)max’  where  W2'  The  shape  of  the  profi1es 
are  similar  in  each  case  but,  as  expected,  the  magnetic  field  which 

reaches  its  maximum  at  closer  distances  to  the  slot,  induced  a 

greater  acceleration  of  the  boundary  layer  fluid  relative  to  that  in 

the  freestream.  The  difference  in  velocity  profiles  is  still  evident 

at  approximately  150  slot  heights  although  its  magnitude  has  negligible 

influence  on  the  film-cooling  effectiveness  (shown  on  Fig.  5.11)  for 

the  same  run  conditions  (curves  1  and  2). 

Figure  5.9  also  shows  how  the  MHD  velocity  boundary-1 ayer 

profiles  become  fuller  with  an  increase  in  the  velocity  ratio  maxima 
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Fig. 5.9  Influence  of  Initial  Condition  for  B(x) on  MHD  Film-Cooling 
Velocity  Profiles 
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Pig. 5.10  Influence  of  Initial  Condition  for  B(x)on  MHD  Film-Cooling  Temperature  Profiles 
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umax^G  father  distances  downstream.  This  is  a  consequence  of  the 
increasing  mass  entrainment  into  the  boundary  layer  flow,  resulting 
from  greater  downstream  electromagnetic  retarding  forces  on  the 
core  flow  and  the  requirement  for  conservation  of  mass  flow  along  the 
channel . 

The  temperature  profiles  resulting  from  the  two  different 

functional  values  for  B  (x)  are  plotted  on  Fig.  5.10.  The  computational 

problems  are  evident  from  the  large  temperature  overshoot  at  y/yP-.5  and 

b 

x/y^=5.  The  profile  at  x/y(_.=5  (for  the  magnetic  field  equal  to  its 
maximum  at  (x/yr)  =1  (dashed  curve))  is  seen  recovering  from  the 
sudden  change  in  magnetic  field.  The  velocity  profiles  in  Fig.  5.9 
(for  the  case  of  magnetic  field  which  is  still  developing,  having 
reached  approximately  50%  of  its  maximum  value  at  x/y^=5)  are  seen  to 
be  well  behaved.  At  twenty  slot  heights,  however,  the  first  profile 
appears  to  have  recovered  (solid  curve)  while  the  second  case  profile, 
that  with  (x/yr)  =10  is  still  showing  signs  of  instability.  It  is 
noted  in  Fig.  5.10  that,  similar  to  the  downstream  convergence  of 
velocity  profiles,  the  temperature  profiles  become  virtually  identical 
at  approximately  150  slot  heights. 

From  these  results,  it  was  decided  that  all  subsequent  com¬ 
putations  would  be  performed  with  a  magnetic  field  B  reaching  its 
maximum  constant  value  at  x/y^=l  since,  among  other  reasons,  the 
magnetic  field  in  actual  MHD  generator  .ducts  would  already  be  fully 
established  at  the  film-cooling  slot.  Any  instabilities  which  were 
encountered  in  the  initial  steps  of  the  computational  marching  pro¬ 
cedure  were  completely  damped  at  twenty  slot  heights  downstream. 
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Fig. 5.11  Influence  of  Initial  Conditioner  B(x)onMHD  Film-Cooling  Effectiveness 


It  should  perhaps  be  noted  that  damping  of  any  initial 
instabilities  could  have  been  accomplished  at  even  smaller  distances 
by  choosing  appropriate  smaller  increments  Ax  in  the  forward  step  of 
the  finite-difference  procedure. 

Figure  5.12  displays  more  intermediate  developing  temperature 
profiles  for  the  same  run  conditions  shown  on  Fig.  5.10;  solid  curves 
represent  B  =B  at  (x/yr)  =10.  It  can  be  seen  that  the  profiles 
converge  at  the  wall  for  downstream  distances  up  to  approximately  30 
slot  heights,  beyond  which  Taw  increases  to  its  recovery  value.  The 
lower  boundary  of  the  mixing  core  emerging  from  the  slot  lip  is 
therefore  seen  to  impinge  on  the  film-cooled  wall  at  x/y^-30  and 
represents  the  point  at  which  the  film-cooling  effectiveness  begins  to 
decrease  from  unity. 
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T/Tg  or  Te/TG 

Fig.  5.12  Developing  Boundary  Layer  Temperature  Profiles  in  MHD  Film-Cooling 
for  B  =  2  at  x/yc  max  =  10  Equilibrium  Case 
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5.4  Thermal  and  Ionizational  Non-Equilibrium  MHD  Film-Cooling 

Solutions  to  Eqs.  (2.4)  and  (2.5),  governing  electron 
conservation  and  temperature,  were  more  difficult  to  obtain.  Because 
of  the  strong  dependence  of  the  plasma  electrical  properties  on  the 
electron  temperature,  T  ,  very  small  forward  steps,  Ax  =  0(10  ),  and 

extended  precision  numerical  computations  as  well  as  an  iteration 
procedure  for  the  wall  electron  concentrations  and  temperatures  were 
requi red . 

The  MHD  influence  on  the  near-slot  velocity  profiles  are 
shown  in  Fig.  5.13,  where  they  are  compared  with  those  of  similar  MHD 
run  conditions  but  for  equilibrium  flow.  The  electromagnetic  body 
forces  for  non-equilibrium  flows  are  seen  to  have  a  similar  influence 
on  the  velocity  profiles  as  was  observed  in  Fig.  5.2  except  that  the 
effect  is  significantly  more  pronounced  for  non-equilibrium  conditions. 
The  larger  j  X  ^  forces  resulting  from  the  relative  higher  current 
densities  (shown  plotted  in  Fig.  5.14)  have  the  effect  of  destroying 
the  velocity  maxima  nearer  to  the  film-cooling  slot  and  cause  the 
profiles  to  become  fuller  at  smaller  downstream  distances  than  was 
the  case  for  equilibrium  MHD  flows,  Fig.  5.13.  The  velocity  profiles 
farther  away  from  the  injection  slot  are  shown  in  Fig.  5.15.  Location 
of  the  velocity  maxima  for  the  non-equilibrium  profiles  in  this 
figure  is  somewhat  lower. 

The  resulting  influence  on  the  near-slot  boundary-1 ayer 
temperature  profiles  can  be  seen  on  Fig.  5.16.  The  temperature  in  the 
boundary  layer  near  the  wall  is  increased  significantly  as  a  result 
of  the  increase  in  higher  temperature  core  flow  becoming  entrained 
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Fig.  5.13  Effect  of  Non-Equilibrium  on  MHD  Film-Cooling  Velocity  Profiles 
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into  the  mixing  region  as  is  evident  from  the  non-equilibrium 
velocity  profiles  shown  previously  in  Fig.  5.13.  It  is  noted,  however, 
that  the  non-equilibrium  effect  on  the  adiabatic  wall  temperature  is 
not  as  pronounced  and  T ^  is  seen  to  be  increased  only  slightly  above 
its  equilibrium  flow  counterpart.  The  electron  temperature  on  the 
other  hand,  has  its  greatest  rise  above  the  main  gas  temperature 
in  the  rear  wall  region  although  any  appreciable  increase  in  heat 
transfer  to  the  wall  due  to  the  higher  electron  temperature  is  not 
evident. 

The  higher  current  densities  in  the  near  wall  region  also 
produce  a  greater  MHD  turbulence  damping  which,  in  turn,  reduce  heat 
transfer  from  the  hotter  far-wall  boundary  layer  region  towards  the 
colder  wall.  The  tangential  injection  of  the  colder  secondary  fluid 
is  seen  to  absorb  any  temperature  increases  near  the  wall  which  would 
be  expected  from  higher  Joule  dissipation  in  the  non-equilibrium  case. 
Farther  downstream  from  the  film-cooling  slot,  however,  some  increase 
in  the  near  wall  boundary  layer  temperature  is  evident  from  Fig.  5.16. 
This  is  due  to  increasing  Joule  dissipation  resulting  from  the  higher 
non-equilibrium  current  densities  as  shown  plotted  on  Fig.  5.14.  The 
magnitude  of  the  maximum  current  overshoot  observed  in  the  present 
solutions  is  considerably  less  than  that  reported  in  Ref.  [7],  where 
the  near  wall  current  densities  were  as  high  as  eight  times  their  free 
stream  value.  It  can  be  expected  that  the  film-cooling  fluid  acts 
to  reduce  the  magnitude  of  the  increase  in  electron  temperature  above 
that  of  the  main  gas,  but  when  compared  to  the  large  current  overshoots 
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Fig.  5.14  Effect  of  Non-Equilibrium  on  Boundary- Layer  Current  Density  Profiles 
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Fig. 5.15  Effect  of  Non-Equilibrium  on  MHD  Film-Cooling  Velocity  Profiles-Far  Slot  Region 
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Fig.  5.17  Influence  of  Non-Equilibrium  onMHD  Film-Cooling  Effectiveness 
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obtained  in  Ref.  [7],  the  significantly  higher  freestream  current 
densities  applied  there  must  also  be  taken  into  consideration. 

The  higher  adiabatic  wall  temperatures  result  in  a  lower 
film-cooling  effectiveness  as  shown  in  Fig.  5.17.  The  effectiveness 
is  lower  than  both  the  equilibrium  MHD  value  and  the  OHD  film-cooling 
effectiveness  in  the  near-slot  region.  At  approximately  60  slot  heights 
downstream,  the  non-equilibrium  effectiveness  becomes  greater  than  its 
OHD  counterpart  although  it  remains  less  than  the  equilibrium  MHD 
value  at  all  downstream  locations.  The  difference  between  equilibrium 
and  non-equilibrium  MHD  effectiveness  also  grows  larger  at  greater 
downstream  distances,  partly  as  a  result  of  increased  non-equilibrium 
current  densities  giving  rise  to  greater  Joule  dissipation. 

The  larger  j  X  B  forces  near  the  wall  can  also  be  seen  in 
Fig.  5.15  to  result  in  a  shift  of  the  maximum  velocity  overshoot 
closer  towards  the  wall  which  in  turn  results  in  a  greater  mixing  of 
the  hotter  free  stream  fluid  with  the  cooling  layer.  Both  this 
effect  and  increasing  Joule  dissipation  near  the  wall  result  in  lower 
MHD  film-cooling  effectiveness  when  compared  with  the  equilibrium  case. 

For  the  cases  computed  (up  to  the  maximum  downstream  station 
indicated  by  Fig.  5.17),  both  equilibrium  and  non-equilibrium  film¬ 
cooling  effectiveness  were  greater  than  their  OHD  counterpart  in  the 
far  downstream  region  which  is  of  most  interest  to  the  designer. 

5.5  Conclusions 

The  fundamental  process  of  tangential -injection  film-cooling 
has  been  investigated  theoretically  for  cooling  the  B-walls  of  both 
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equilibrium  and  non-equilibrium  MHD  generators.  The  general  finite- 
difference  method  of  Patankar  and  Spalding  for  solving  heat  and  mass 
transfer  boundary  layer  flows  has  been  modified  in  a  universal  way  to 
allow  MHD  heat  and  mass  transfer  flow  analysis.  The  free  molecular 
sheath  analysis  required  for  non-equilibrium  conditions,  however, 
requires  individual  analysis  for  the  flow  conditions  under  consider¬ 
ation  . 

The  overall  influence  of  electromagnetic  body  forces  on  the 
adiabatic  film-cooling  effectiveness  results  in  a  significant  differ¬ 
ence  over  its  OHD  counterpart.  As  the  magnetic  field  induction,  B, 
is  increased  from  zero  to  moderate  values,  an  increase  in  n  is  observed 
and  this  increase  grows  larger  with  increasing  distance  downstream 
from  the  coolant  injection  slot.  As  an  example  for  the  flow  conditions 
shown  in  Fig.  5.3  an  increase  of  23%  is  observed  at  x/y^=100  and 
86%  at  x/y^=500  when  B  is  changed  from  0  to  1 .  As  the  magnetic  field 
induction  is  increased  further  to  B=2  the  film-cooling  effectiveness 
is  not  significantly  altered  farther  downstream.  Nearer  the  slot, 
however,  q  is  decreased  and  for  x/y^  <^60  the  decrease  results  in  a 
lower  film-cooling  effectiveness  when  compared  with  the  OHD  case.  As 
B  is  further  increased  to  extremely  large  values,  the  film-cooling 
effectiveness  continues  to  decrease  in  the  near-slot  region  whereas 
farther  downstream  boundary  layer  separation  or  choking  of  the  channel 
flow  will  occur.  Thus,  there  appears  an  optimum  value  of  magnetic 
field  intensity  which  will  result  in  maximum  film-cooling  effectiveness 
for  any  one  particular  set  of  flow  conditions.  For  the  run  conditions 
of  Fig.  5.3,  this  optimum  magnetic  field  occurs  for  values  of  B 


between  1  and  2. 
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Allowance  for  electromagnetic  field  damping  of  turbulent 
fluctuations  has  an  appreciable  influence  on  the  MHD  film-cooling 
effectiveness.  With  a  particular  value  of  B=2,  for  example  and  the 
flow  conditions  indicated  on  Fig.  5.3,  MHD  damping  of  turbulent  fluc¬ 
tuations  results  in  an  average  15%  increase  in  the  film-cooling  effec¬ 
tiveness.  Allowance  for  MHD  turbulence  damping  should  therefore  be 
included  in  any  MHD  film-cooling  predictions. 

The  influence  of  initial  flow  parameters  and  injection  slot 
geometry  on  the  film-cooling  effectiveness  is  not  the  same  for  MHD 
flows  when  compared  with  OHD  flows.  Because  the  velocity  profiles  are 
significantly  altered  in  MHD  boundary  layers,  the  optimum  film-cooling 
effectiveness  is  not  necessarily  obtained  with  a  cool ant-to-primary 
flow  velocity  ratio  of  unity  as  has  been  shown  the  case  in  OHD  film¬ 
cooling  processes.  The  optimum  velocity  ratio  Uq/Uq  to  achieve  maximum 
effectiveness  at  any  given  magnetic  field  intensity  has  not  been 
investigated  in  the  present  work.  This  could,  however,  be  determined 
if  many  more  runs  were  performed  with  the  present  numerical  method. 

The  MHD  film-cooling  effectiveness  is  influenced  to  a  lesser 
degree  than  its  OHD  counterpart  for  similar  changes  in  the  initial 
velocity  and  temperature  ratios ,  Uc/Uqo  and  T^/T^.  .  The  MHD  film¬ 
cooling  effectiveness  does,  however,  depend  on  the  magnitude  of  these 
initial  profiles.  For  a  constant  magnetic  field  induction,  increases 
in  effectiveness  are  relatively  larger  at  lower  values  of  initial 
velocities  and  temperatures .  A  definite  range  of  magnitudes  for  initial 
velocities  and  temperatures  is  obtained  which  result  in  negligible 
further  changes  in  effectiveness  at  any  given  magnetic  field  intensity. 
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These  observations  were  made  by  keeping  the  ratios  uq/ Uq0  and  T^/Tg 
constant  at  0.5.  The  influence  of  variable  initial  ratios  at  constant 
magnetic  field  has  not  been  investigated  here  although,  again,  such 
sensitivity  analysis  could  be  performed  with  additional  runs  of  the 
present  computational  procedure. 

The  sensitivity  of  the  MHD  film-cooling  effectiveness  to 
initial  flow  conditions  is  no  more  pronounced  than  what  has  been 
observed  in  OHD  film-cooling  studies.  The  choice  of  initial  condition 
for  the  magnetic  field  induction,  however,  influences  the  downstream 
solution  for  n  in  a  more  significant  way.  Although  the  percentage 
difference  in  magnitude  observed  in  the  effectiveness  is  relatively 
small,  this  difference  is  sustained  for  longer  distances  downstream 
from  the  injection  slot.  It  is  therefore  more  important  to  perscribe 
initial  magnetic  field  variations  as  accurately  as  possible  than  is 
necessary  for  other  initial  flow  conditions  in  either  MHD  or  OHD  film- 
cool  ing  predictions . 

Thermal  non-equilibrium  and  finite-rate  ionization  have  a 
significant  influence  on  film-cooling  effectiveness.  The  major  cause 
for  this  is  the  alterations  of  the  boundary  layer  velocity  profiles 
resulting  from  larger  J  X  ^  forces.  These  in  turn  are  the  result  of 
increased  near-wall  current  densities,  je  ,  due  to  the  strong  depen¬ 
dence  of  plasma  electrical  properties  on  electron  temperature  and 
number  density. 

In  the  near-slot  region  the  coolant  fluid  precludes  any 


significant  thermal  non-equilibrium  and  temperature  rise  which  would 
be  attributed  to  Joule  dissipation.  There  is  a  bulk  temperature  rise 
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in  the  mixing  region  as  a  result  of  the  flow  velocity  distribution 
although  the  adiabatic  wall  temperature,  and  hence  the  film-cooling 
effectiveness,  is  changed  only  slightly  from  that  in  equilibrium  MHD 
flows.  Farther  downstream,  however,  an  increasing  current  density 
over-shoot  develops  near  the  wall  which  influences  the  film-cooling 
effectiveness  in  several  ways.  First,  a  larger  velocity  maxima  is 
produced  in  the  boundary  layer  and  the  resulting  wall-jet  like  flow 
serves  to  increase  the  film-cooling  effectiveness  in  a  manner  similar 
to  the  one  observed  in  the  equilibrium  case.  The  shifting  of  this 
velocity  maxima  closer  to  the  wall,  on  the  other  hand,  brings  more  of 
the  hotter  core-flow  gases  in  closer  contact  with  the  coolant-gas  layer 
producing  a  more  rapid  temperature  rise  of  the  coolant  gas.  Secondly, 
a  larger  current  density  produced  greater  MHD  damping  of  turbulent 
fluctuations  such  that  the  rate  of  this  temperature  rise  is  tempered 
to  some  degree. 

Finally,  the  increasing  current  densities  at  larger  down¬ 
stream  distances  resulting  from  greater  non-equilibrium  conditions, 
produce  a  slight  increase  in  the  near-wall  temperature  due  to  Joule 
dissipation.  This  increase  is  noticeable  but  relatively  small  for  the 
flow  conditions  assumed.  The  net  influence  of  non-equilibrium 
therefore  is  to  produce  a  film-cooling  effectiveness  that  is  reduced 
from  its  equilibrium  counterpart,  but  remains  greater  than  that  for 
an  equivalent  OHD  flow. 

For  equal  wall  thermal  protection,  therefore,  more  film¬ 
cooling  slots  are  required  for  non-equilibrium  MHD  generators  than  for 
equilibrium,  combustion  driven  generators.  As  a  matter  of  relative 
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comparison,  if  the  film-cooling  effectiveness  for  a  channel  wall  were 
permitted  by  design  to  reduce  to  50%  of  its  initial  value,  then  for  a 
particular  set  of  flow  conditions,  the  MHD  equilibrium  channel  should 
have  tangential -i njection  spaced  at  250  slot  heights,  the  MHD  non¬ 
equilibrium  channel  at  190  slot  heights  and  the  OHD  channel  at  120 
slot  heights. 
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Abstract 

The  Karman-Pohl hausen  technique  is  applied  to  obtain  a  solu¬ 
tion  of  MHD  boundary  layer  flow  over  an  insulated  semi -i nf i ni te  flat 
plate  in  the  presence  of  a  uniform  transverse  magnetic  field  and  uni¬ 
form  suction  applied  at  the  plate  surface.  The  fluid  assumed  to  be 
incompressible  and  have  a  constant  coefficient  of  electrical  conduc¬ 
tivity.  The  usual  MHD  assumptions  are  adopted  among  which  the  magnetic 
Reynolds  number  is  considered  to  be  small  enough  such  that  induced 
current  effects  become  negligible.  Two  cases  are  considered:  the 
first  is  usually  referred  to  as  magnetic  field  fixed  relative  to  the 
fluid  far  from  the  plate  and  the  second  designated  as  magnetid  field 
fixed  relative  to  the  plate.  For  the  first  case  an  asymptotic  -  MHD 
suction  profile  is  obtained  and  the  advent  of  the  fully  developed 
boundary  layer  is  illustrated  for  combined  Hartmann  numbers  and  suction 
coefficients.  Developing  velocity  profiles  and  average  skin  friction 
coefficients  are  also  derived.  For  the  second  case,  boundary  layer 
separation  is  predicted  for  an  impermeable  wall  at  a  Stuart  number  of 
unity  in  spite  of  a  nonexistent  pressure  gradient.  Unfortunately,  the 
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integral  method  with  a  fourth  order  polynomial  velocity  profile  was 
not  able  to  predict  the  points  of  separation  for  suction  and  magneti 
field  applied  simultaneously. 
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Nomencl ature 


B 

c, 


L 

M 

N 

Nx 

Re, 

Re) 

Ue 

U. 

u 


V 

w 

x,y 

6 


6 


1 


6 


2 


v 


P 

a 


x 


w 


=  magnetic  field 

=  local  skin  friction  coefficient  (=x  /A  pU 
=  average  skin  friction  coefficient 
=  suction  coefficient  (=|v  l/U  ) 

=  reference  plate  length 
=  Hartmann  number  (=  BL/a/y^) 

2 

=  interaction  parameter  (=  aB  L/pUj 

=  Stuart  number  or  local  interaction  parameter  (=aB  x/pUj 
=  Reynolds  number  (=UcoL/v) 

=  local  Reynolds  number  (=Uoox/v) 

=  free  stream  velocity 
=  uniform  upstream  velocity 
=  x-component  of  velocity 
=  uniform  suction  velocity 

=  coordinates  parallel  and  normal  to  plate,  respectively 
=  boundary  layer  thickness 
=  displacement  thickness 
=  momentum  thickness 
=  fluid  dynamic  viscosity 
=  fluid  kinematic  viscosity 
=  fluid  mass  density 
=  fluid  electrically  conductivity 
=  wal 1  shear  stress 


Subscripts 

B  =  magnetic  field 

e  =  free-stream  conditions 

r  =  reference  based  on  plate  length  L 

s  =  separation 

w  =  evaluated  at  the  wall 

0  =  zero  suction 

00  =  uniform  upstream  conditions 
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Introduction 

Magnetohydrodynamic  (MHD)  flows  over  a  porous  surface  are  en¬ 
countered  in  studying  vortex  type  MHD  generators^  or  when  considering 
the  possibility  of  applying  transpiration  cooling  to  the  channel  walls 
of  an  MHD  converter.  In  such  cases  it  becomes  necessary  to  solve  the 
flow  equations  in  the  presence  of  permeable  wall  boundary  conditions. 
Hydromagnetic  flows  with  porous  walls  have  interested  previous  authors 
among  which  are  E.R.  D'Sa  who  obtained  a  series  solution  for  the  hydro- 
magnetic  wall  jet  including  the  effects  of  suction  and  blowing  at  the 
wall  and  R.K.  Jain  and  K.N.  Mehta  considered  the  problem  of  laminar 
hydromagnetic  flow  in  an  annulus  with  porous  walls.  In  general  the 
problem  is  complicated  but  simplifying  assumptions  are  possible  that 
permit  solutions  from  which  first  order  effects  can  be  inferred. 

The  Karman-Pohl hausen  method  has  already  been  found  useful 
in  solving  the  MHD  boundary  layer  equations  in  the  absence  of  suction 

4 

or  blowing.  Maciulaitis  and  Loeffler  used  a  parabolic  velocity  profile 

5 

in  their  investigation  of  MHD  channel  entrance  flows,  W.C.  Moffatt 

chose  a  profile  that  was  a  function  of  the  Hartmann  number  in  his 

study  of  hydromagnetic  channel  entrance  flows.  Heiser  and  Bornhorst6 

proposed  a  modified  fourth  degree  polynomial  velocity  profile  to  include 

adverse  free-stream  pressure  gradients  for  flows  of  incompressible, 

conducting  fluids  over  semi -infinite  boundaries.  The  modification 

relinquishes  the  usual  ordinary  hydrodynamic  (OHD)  condition  that 
2  2 

3  u/3y  vanish  in  the  free  stream  and  replaces  this  by  a  condition  on 
the  third  derivative  at  the  wall.  Later,  Ramamoorthy^  resolves  the 

o 

problems  treated  by  Rossow  and  others  using  the  Karman-Pohl hausen 
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method  with  the  modified  velocity  profile  and  obtains  excellent  agree¬ 
ment  as  well  as  improvements  over  the  solutions  obtained  earlier.  In 
the  present  investigation  the  modified  fourth  degree  polynomial  profile 
is  not  necessitated  since  an  adverse  pressure  gradient  is  not  considered. 
For  the  case  referred  to  by  Rossow  as  "magnetic  field  fixed  relative 

to  the  fluid",  an  exponential  profile  similar  to  that  prescribed  by 
9 

H.G.  Lew  for  OHD  flat  plate  flow  with  uniform  suction  will  be  employed. 
For  the  second  case  referred  to  as  "magnetic  field  fixed  relative  to 
the  plate",  the  usual  fourth  degree  polynomial  profile  is  adopted  since 
in  general  this  yields  slightly  better  results  than  the  exponential 
assumption. 

Governing  Equations  and  Solutions 

Our  mathematical  model  consists  of  a  boundary  layer  flow  over 
an  electrically  insulated  flat  plate  in  the  presence  of  a  constant 
transverse  magnetic  field  B  which  may  be  stationary  or  moving  with 
constant  velocity  Ug  relative  to  the  plate.  The  fluid  is  assumed  to 
be  incompressible  with  constant  electrical  conductivity  a,  density  p 
and  viscosity  v.  We  further  assume  that  no  electric  field  is  applied 
and  that  the  uniform  upstream  velocity  is  undisturbed;  that  is,  the 
magnetic  field  is  assumed  to  vanish  abruptly  upstream  from  the  leading 
edge  of  the  plate. 

Application  of  the  integral  method  results  in  the  following 
general  differential  equation  which  is  comparable  to  that  given  in 
References  4,  5  and  6  in  the  absence  of  wall  suction: 
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dx 
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The  first  case  of  interest  arises  when  UD  =  U  =  U  are  constant. 

B  e  °° 

g 

Following  Lew  we  introduce  new  variables  X  and  E,  defined  by 
A(S)  =  -  ReL  CQ  6r(£)/2  and  E,  =  Re[_2  CQ2  x  =  CQ2  Rex  with  =  6/L, 
x  =  x/LRe^,  y  =  y/6  and  u  =  u/U^  we  obtain  the  following  dimension¬ 
less  equation: 
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where  the  combined  magnetic  and  suction  effects  are  related  by  a  singl 
parameter  I  defined  by  I  =  N/Re^Cg  =  (M/Re^Cg)  .  Even  for  an  imperme 
able  wall  the  boundary  layer  behaves  asymptotically  as  x  ->  °°.  This 
asymptotic  solution,  as  given  by  Gurevich  et  al^,  is 


u  =  1 


(3) 


where  yr  =  6 =  y/L.  To  solve  equation  (2)  we  will  want  to  choose 
an  exponential  profile  of  a  particular  form  and  to  illustrate  the 
justification  for  this  choice  we  first  resolve  the  momentum  equation 
for  zero  suction  and  compare  the  results  with  several  existing  solu¬ 
tions.  The  non-dimensional  equation  to  be  solved  is  now 


d£ 


[X 


0 


0 


u  ( 1  -ID  dyj  +  (^-)  A^  (l-u)dy  = 


0 


1  (^) 


o 


4A0  3y 


(4) 


w 


2— 

with  Aq  =  Re^6r/2  and  ^  =  Re^  x.  Considering  the  asymptotic  solution 
(3)  we  are  led  to  assume  a  velocity  distribution  of  the  form 
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-My 

u  =  1  -  e  (l-MyK(iT)) 


(5) 


1  2 

where  K(x)  =  ^  ( <S  (3T)  -  1)  as  determined  from  the  boundary  condition 
(9  u/9y  )w  =  -  M  6  .  Expression  (5)  satisfies  the  remaining  boundary 
conditions  u(0)  =  U'(°°)  =  u''(oo)  =  u(«)  -1=0  and  we  note  that  when 
<5r  =  1  ,  K  =  0 ,  y  =  y^  and  hence  the  required  asymptotic  solution  (3) 
is  recovered.  Substitution  of  (5)  into  equation  (4)  and  defining 
=  resu^ts  in  the  differential  equation 


d^n  m  2 


84)q4  -  10c|>02  +  1 

L  ^q(  ~44)q4-1  6cfQ3-l  24>02+l  1 ) 


(6) 


The  solution  of  (6)  subject  to  =  0  at  =  0  is  easily  found  to  be 


^Re^  ^0  Nx  "  8  ^0  '  ^ 


0 


iVln{ 


t  C  C  I 

(4>0+ai )  (ar40)  ( <t)o+a2 ^ 

^5  C 
a1  0  a2L 6 


(7) 


with  the  constants  given  by 


17a13-20a12-n^-  a-,+2 
2a i  (a-j  2-a22) 

-17a23+20a22+ll^-  a2-2 
2 ^ C ^ "j  ”a2  ) 


1  7a-,  3+20a-,  2-l  1  t>-  3^2 
2^1  ( a i  ~^2  ^ 

-17a23-20a22+ll2-  a2+2 
2a2(a-j2-a22) 
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Cr  = 


al  = 


34a-j  -23 
2(a12-a22) 

/  5+/1 7v 1  /2 
1  8  1 


> 


-34a2  +23 
2  27 

2  ( a -j  -a 2  ) 


=  ( 


5-/17 

8 


1/2 


) 


Using  solution  (7)  the  local  skin  friction  coefficient  and 

displacement  thickness  can  now  be  plotted  as  a  function  of  Stuart 

number  N  given  that 
x  J 


Re^ 

lA 


40  ^0  +  2) 


(8) 


and 


61 M  ^  2 

=  24>0(f  -  2<J>g) 


(9) 


Figure  1  shows  fairly  good  agreement  with  the  other  results  shown  al¬ 
though  a  somewhat  higher  skin  friction  coefficient  is  obtained.* * 

Returning  now  to  the  solution  of  equation  (2)  with  suction 
present,  an  asymptotic  -  MHD  suction  profile  is  easily  shown  to  be 

u  =  1  -  e'An  (10) 

where  A  =  j  +  /(^  +  I)  and  n  =  y | v  | /v .  An  obvious  choice  for  u 
now  is  expression  (5)  with  A  simply  replacing  the  Hartmann  number  M. 
Hence  we  assume 

*  7 

It  appears  that  Ramamoorthy  incorrectly  plotted  c|M  instead  of 

*8 

c^Re^/M.  Expressions  used  to  plot  Rossow's  results  for  Nx  £  .2  are 

c^Re./M  =  ( .664+2. 29N  -2. 768N  2)//N  and  6-.M/L  =  (1 .73-2.21N  )/N 
T  L  XXXI  XX 


. 


SKIN  FRICTION  COEFFICIENT 
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isplacementThickness  and  Skin-Friction  Coefficient 
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Fig.  2  Velocity  Profiles  at  Various 
Locations  and  Along  Plate  Length 


-Ay 

u  =  1  -  e  (l-AyK(x)) 


(ID 


The  parameter  K(3T)  is  then  given  by 


K(x) 


[A(A+2A)-4IA2] 

2A(A+A) 


(12) 


obtained  from  the  condition  that  at  the  plate  surface  -  Re,  Cn  6 

L  (J  r 

_  2 _  _ 2  2 

(du/3y)w  =  (9  u/3y  )w  +  N  ReL  6^.  The  remaining  boundary  conditions 
are  again  automatically  satisfied  by  expression  (11).  Furthermore, 
the  asymptotic  solution  is  also  recovered  with  K  =  0,  which  occurs 
when  A  =  (A-A/1+4I ) / 4 I  =  -  ^  >  since  with  | A |  =  i,  y  =  n  in  view  of 
n  =  2 | A | y .  Insertion  of  (11)  into  the  Kdrmdn  equation  (2)  results 
in  the  ordinary  differential  equation 


dA  _  2A(A+A)2[(2AA+A2)2  +  81  A2(2I A2-2AA-A2) ] _ 

dC  12A2A4+36A3A3+37A4A2+11A5A-8IA3(6IA3+10AIA2+A2A+3A3) 

g 

which  agrees  with  the  non-magnetic  case  of  Lew  for  I  =  0  (A=l) 

2 

specialized  for  zero  Mach  number  and  with  the  coefficient  of  A  corrected 
to  37.  Equation  (13)  is  not  easily  integrated  in  closed  form  for  I  f  0 
and  a  numerical  solution  was  obtained  subject  to  A  =  0  when  £  =  0. 

The  developing  velocity  profiles  for  several  values  of  the 
magnetic  suction  parameter  I  are  shown  in  Fig.  2.  The  displacement 
thickness,  defined  by 


6 


1 


A  r  4AA+3A2-4IA2  i 
pZ  1  A+A  * 


v 


(14) 


g.  3  Displacement  thickness  along 
plate  length,  fixed  to  fluid 
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Fig.  4  Position  of  fully  developed  MHD 
boundary  layer  with  uniform  suction, 
F  fixed  to  fluid 


139 


is  plotted  in  Fig.  3.  It  is  evident  that  the  displacement  thickness  is 
reduced  and  the  fully  developed  boundary  layer  is  approached  at  smaller 
£  as  I  increases.  Expression  (14)  is  easily  verified  to  approach  its 
asymptotic  value  of  1/A  when  X  =  -  1/2.  Defining  the  fully  developed 
condition  to  occur  at  99%  of  this  value  it  is  possible  to  illustrate 
the  combined  effect  of  suction  and  magnetic  field  on  the  attainment 
of  the  asymptotic  state.  Figure  4  indicates  the  position  on  the  plate 
corresponding  to  fully  developed  conditions  as  a  function  of  M/Re^ 
(ratio  of  ponderomotive  force  to  inertia  force)  with  Cq  as  a  parameter. 
For  sufficiently  weak  magnetic  fields  (or  low  electrical  conductivity) 
the  magnitude  of  suction  primarily  determines  the  asymptotic  position. 
Conversely,  for  large  values  of  M/Re^  the  behaviour  corresponds  more 
as  if  the  plate  were  impermeable.  Gurevich  et  al^  determine  from  a 
series  solution  without  suction  that  at  a  Stuart  number  of  approxi¬ 
mately  2  or  larger,  the  boundary  layer  behaviour  differs  little  from 
its  asymptotic  state.  From  the  present  analysis  a  corresponding  value 
of  N  =  1.42  has  been  obtained. 

A 

It  is  also  desired  to  know  the  influence  of  magnetic  field 
and  suction  on  the  average  skin  friction  coefficient  defined  as 


C.  = 


j  pU  2x 

2  00 


/  3u  \ 

yf  V 


dx  = 


w 


.  ^ 

K 


}(^)  dc 

3*w 


(15) 


Figure  5  shows  the  behaviour  of  Cp  with  local  Reynolds  number.  The 
asymptotic  solution  of  2ACq  is  indicated  as  well  and  the  approximate 
Reynolds  number  at  which  it  is  approached  can  be  determined  from  Fig.  4. 


103M/ReL  = 
31.6  |I  =  10) 
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Fig.  5  Average  Skin- Friction  Coefficient,  B  fixed  to  fluid 
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O 

Also  shown  are  some  zero  suction  results  obtained  by  Rossow. 

The  second  case  of  interest  occurs  when  Ug  =  0.  The  resul- 

O  11 

ting  differential  equation  is  that  considered  by  V.J.  Rossow  *  and 
2 

E.R.  D'Sa  .  Assuming  that  no  magnetic  interaction  occurs  with  U  ,  the 

p 

free  stream  velocity  is  now  given  by  Ue  =  -  aB  x/p  or  Ue/Uoo  =1  -  IE,. 

The  dimensionless  momentum  integral  equation  is  now 


1 

2(1  -I?)  ^j={A 

o 


fl 

u*(l-u*)dy]  -  41 X 

0 


u*(l -u*)dy 


+  (— )  /2A 


9y 


w 


(16) 


where  u*  =  u/Ue.  Since  in  this  case  there  is  no  a  priori  reason  to 
assume  a  velocity  profile  of  a  particular  form,  the  usual  polynomial 
distribution  is  adopted-viz., 


u* 


_  _9  — 3 

ay  +  by  +  cy  +  dy 


(17) 


with  coefficients  given  by  a  =  6/(A+3),  b  =  6A/(A+3),  c  =  -  2(4A+3)/(A+3) 
and  d  =  3(X+1 )/ (X+3) .  Substitution  of  (17)  into  equation  (16)  results 
in  the  equation 


dX  _  1 

d^  '  2(1-10 


35(X+3)2(X2+3X+3)  +  4IX2(X+3)(4X2+26X+37) 

- j - 3 - p -  } 

4X4  +  36XJ  +  119X^  +  1 1 IX 


which  once  again  can  be  made  to  agree  with  the  incompressible  version 

Q 

of  Lew's  equation  (19)  when  1=0. 

The  displacement  thickness  and  local  shear  stress  coefficient 


are  given  by 


(19) 


w  A(4A+9) 
1  v  "  5( A+3) 


and 

6Cq(1-U) 
f  "  A(A+3) 


(20) 


An  unfortunate  shortcoming  of  assumption  (17)  is  that  the  asymptotic 
state  for  the  non-magnetic  case  cannot  be  attained.  As  x  + 

5 1 1  Vw I / v  ^  anc*  ecluatl'on  (19)  shows  that  A  for  this  value  is  imaginary 
Lew  has  shown  that  the  useful  range  of  values  for  A  is  only  0  >_  A  >  -  1 
which  also  turns  out  to  be  the  case  for  I  f  0. 

Again  we  resolve  the  purely  magnetic  case  using  (17)  without 
the  modification  proposed  in  Ref.  6,  and  compare  the  results  with  an 
existing  solution.  The  non-dimensional  equation  to  be  solved  for  this 
case  is 


1 

n  -  [Ao  f 

o 


1 


2<r^ 


u*(l -u*)dy  = 


i  (M) 


4A0  3y 


(21) 


w 


and  the  coefficients  in  (17)  take  the  values  of  a  =  2,  b  =  0,  c  =  -  2 
and  d  =  1.  The  solution  for  the  Stuart  number*  is  now 


N 


x 


[1 


148 

315 


(M<t>0)2]1/4 


(22) 


The  Stuart  number  N  is  identical  with  the  parameter  mx  defined  by 
Rossow8.ll.  x 
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and  the  expression  for  the  average  skin  friction  becomes 


C 


F 


M  148 
ReL  315 


*0 

/(I  -  (MV2) 


(23) 


A  comparison  of  Cp  obtained  from  the  polynomial  distribution  with  that 

o  o 

obtained  by  Rossow  with  M/Re^  =  10"  is  shown  in  Table  1 .* 


Nx 

Rossow** 

Eqn.  (23) 

0.01 

0.0132 

0.0133 

0.10 

0.0038 

0.0040 

0.15 

0.0029 

0.0032 

0.20 

0.0024 

J 

0.0026 

Table  1  Comparison  of  Average  Skin  Friction 
Coefficient  Cp,  M/Re^  =  10"^ 

Equation  (18)  was  again  integrated  numerically  satisfying 
the  condition  X  =  0  at  £  =  0.  The  local  skin  friction  from  equation 
(20)  is  plotted  in  Fig.  6  at  a  particular  location  on  the  plate  given 

4 

by  Re  =  10  .  The  restricted  range  of  X  is  seen  to  be  a  limitation  on 

A 

the  magnetic  parameter  M/Re^  for  which  a  solution  is  obtained  when 
suction  is  present  and  a  continuation  up, to  the  point  of  separation 
(c.p  =  0)  was,  therefore,  not  possible.  Nevertheless,  the  combined 

*  1. 328-1. 19N  "  1 .328-1 . 1 9Ny 

Rossow's  result  is  Cr  =  - =  n—  {  — ■ — - -  },  Nv  <  .2 

*■  /  Re  KeL  /  Nv  x  “ 

A  A 
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Fig.  6  Local  skin-friction  variation  at 
a  particular  location  given  by  Rex  =  10  , 
B"  fixed  to  plate 


effect  of  suction  and  magnetic  field  on  the  local  shear  stress  before 
separation,  is  indicated.  When  Cg  =  0,  however,  the  expression  for 
the  local  shear  stress 


Re 


L  


2d-Nx) 


f  M 


r31  5  £ _ 1 


148  (1-NX)4 


-  I]) 


1/2 


(24) 


indicates  that  separation  occurs  when  the  Stuart  number  N  =  1,  which 

X 

_2 

in  Fig.  6*  occurs  when  M/Re^  =10  .  If  a  dimensionless  point  of 

A 

separating  is  defined  by  xs  =(Uoo/v)xs  we  can  plot  the  relationship 
between  k  and  x$  where  k  is  a  constant  in  the  free-stream  potential 

flow  of  the  form  U  /U  =1  - ( kU  /v)  x.  Figure  7  shows  x„  as  a  function 

e  oo  oo'  3  s 

A 

of  k  and  is  compared  with  x$  when  separation  is  due  to  an  adverse 
pressure  gradient.  In  the  latter  case  the  well  known  relationship  in 
Reference  12  is  x$  =  0.125/k  whereas  from  the  present  analysis  we  ob¬ 
tain  xs  =  1/k.  The  reason  for  the  downstream  displacement  of  the 
separation  point  for  the  magnetic  case  is  a  result  of  the  fact  that  a 
hydrodynamic  pressure  gradient  is  transmitted  to  the  plate  surface 
whereas  the  magnetic  body  force  vanishes  there  for  this  case. 

It  would  be  of  interest  to  include  suction  or  blowing  in  the 
analysis  of  channel  entrance  region  problems  considered  in  References 

4  and  5.  If  a  polynomial  profile  were  assumed  the  modification  pro- 

3—3  2 

posed  by  Heiser  and  Bornhorst  would  require  (9  u*/3y  =  NRe,  <5  (9u*/9y) 

W  Li 

-  CQReL6r(32u*/3y2)w. 


* 

Rossow's  result  shown  for  comparison,  was  terminated  at  M/Re^  ^  4.5  x 
10-3  since  it  ceases  to  be  valid  beyond  N  =  0.2. 

X 


PRESSURE  GRADIENT  WITHOUT 
MAGNETIC  FIELD,  REF.12 

PRESSURE  GRADIENT  WITH 
MODIFIED  POHLHAUSEN 
METHOD,  REF.  6 

ZERO  PRESSURE  GRADIENT 
WITH  MAGNETIC  FIELD, 
PRESENT  WORK 
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Fig.  7  Dimensionless  point  of  separation  for 
uniformly  decelerated  flows  of  the 
form  U  /U  =  1  -  kx 

e  oo 
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